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Preface 


Fourier Analysis has now become a widely used technique in many 
branches of science and engineering. This book has been written 
mainly with the needs of students of these subjects in mind. Thus full 
mathematical rigour has not been attempted. On the other hand an 
attempt has been made to avoid presenting the subject in such a way 
that it becomes a mere 4 tool ’. An elementary proof of the fundamental 
theorems has been given to provide some insight into the more 
theoretical aspects of the subject and to provide a better basis for the 
understanding of some of the phenomena which occur later, such as 
the Gibbs’ phenomenon. The first two chapters deal with Fourier 
Series. These are followed by two on Fourier Integrals and finally two 
on applications. It is common to think of Fourier Analysis in terms 
of the analysis of time varying functions. The final chapter shows that 
this point of view is unnecessarily restricted and that the theory is of 
much wider applicability. 

In compiling a work of this nature the work of many other authors 
has of course been drawn upon and in particular H. S. Carslaw, 
Theory of Fourier Series and Integrals should be mentioned. Other 
works are acknowledged in the text. The author’s views on the subject 
were influenced to a considerable extent also by a series of lectures 
given by Mr J. A. Ratcliffe at the Cavendish Laboratory in 1948. 
Finally, the author would like to thank Professor R. D. Klein for 
reading the manuscript and making many helpful criticisms and 
suggestions. 
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CHAPTER I 


Fourier Series 


1.1 Introduction 

Our studies in acoustics introduce us to the idea of harmonics. A 
string fixed at each end can vibrate at a number of different frequencies. 
The lowest of these, /, is called the fundamental frequency. The 
others will have values 2/, 3/, etc., and are called harmonics. These 
frequencies can be excited simultaneously and the resulting vibration 
has a complex waveform. This will still be periodic and the period T 
is equal to (1 If). Thus we see that the summation of a number of 
frequencies harmonically related gives rise to a complex periodic 
waveform. Conversely a complex periodic waveform can be analysed 
into a number of sinusoidal variations which are harmonically related. 
It is the second process we shall be concerned with in this chapter. 
The plot of amplitude against frequency which in this case consists of 
a number of discrete lines at frequencies /, 2 f 9 ... is called the spectrum 
of the waveform. 

1.2 The Fourier Series 

Suppose we have a function /(/) which is periodic and suppose for 
simplicity that the period is 2n f i.e./(f+27r) =/(f). Then/(f) can be 
represented as a series 

/(/) = ia 0 +aiCost+a2COs2t+... + b l sint+b 2 sin2t+... 

00 

= i a 0 + 2 (fl rt cosrt/+& n sin/ 2 /) 

n= 1 

This is the Fourier series representation of the function /(f). As 
might be expected certain restrictions must be placed on /(f) for the 
expression to be valid. The conditions which must be satisfied, known 
as Dirichlet’s conditions, are met by practically all of the functions 
which are of interest in the physical sciences. The integral J%|/(f)| dt 
must be finite and /(f) must be piecewise continuous and piecewise 

11 
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monotonic. If/(f) is piecewise continuous then its discontinuities are 
restricted to a finite number of jump discontinuities in any finite 
interval of f, a jump discontinuity being one in which the change in 
/(f) as we pass through the discontinuity is of finite magnitude. If/(f) 
is piecewise monotonic then it is possible to divide any finite range of 
f into a finite number of intervals in such a way that within each inter¬ 
val /(f) is either monotonically increasing or monotonically decreas¬ 
ing (i.e./(f) has no maxima or minima within the interval). We shall 
rarely have occasion to refer to these conditions in the present work. 

The expression above for/(f) represents the function/(f) for all f. 
Evidently if we had a function/(/) which was not periodic but which 
was defined only over the finite range of f, - tt < t < tt then the 
expression would represent/(f) for the values of f for which/(f) is 
defined, i.e. — tt ^ f ^ tt. 

More generally if the function /(f) is periodic and has period T 
then the expansion becomes 


fit) = K+2 cos y sin y 


A non-periodic function defined over a finite range - £7^ f ^ \T 
is also represented by this expression. Outside the range in which/(f) 
is defined, i.e. |f| > \T the expression represents a periodic extension 
of the function /(f). 


1.3 Properties of Sine and Cosine 

We now list some important integral properties of sine and cosine 
functions which we shall need and which can readily be verified using 
elementary calculus. If/? and q are integers (other than zero) then we 
have 

j^cosptdt = 0 
j^sinptdt = 0 
j”^ cos 2 ptdt = TT 
J^sin 2 ptdt = tt 
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J ^oospt.cosqtdt = 0 , p ^ q 
j”^sinpt.s'mqtdt = 0 , p ^ q 
J _^cospt.smqtdt = 0 

The last three of these relations show that sines and cosines form an 
orthogonal set of functions. (They include the first two relations which 
are just special cases.) 


1.4 Evaluation of the Coefficients 

Using the relations given in the previous section we can evaluate the 
coefficients a m b n . We have 

/(f) = + o j cos f+ 2 cos 2f+... + £> j sin f+ 62 sin 2f +... 

Hence integrating over the period — tt to tt and using the relations 
of section 1.3 we have 


/: 


fit) dt = la 0 .2n 


ao 


-if 

KJ-T 


mdt 


Similarly multiplying through by cos nt and integrating we have 
/(f)cos ntdt = a n TT a m = - j /(f)cos ntdt 

j—v I* J — TT 

Multiplying through by sin/jf and integrating we have 

I /(f)sin//f dt = b n TT .*. b n = - f 

J-fT TTj-TT 

For the more general case 

2 r 

•• ao= rJ_ 


/(f) sin wf dt 


/(f) dt 


-Til 
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rri2 

J-772 


/WC 0 SyW/<// 


/. a 


n 


and 



„ v . 2ir _ 
f(t)sm — ntdt 


/: 
is: 
s: 


T/2 2 It 

o n cos 2 —ntdt = |aj 

-772 * 

2 r 772 v ^ 7 T 

f{t) cos—ntdt 
-r/2 


T/2 27T 

sin 2 —*/<// = 
-r/2 1 


. b _ 2 r 


2 f r/2 2 tt 

/(Osin— «f<// 
r/2 ** 


These expressions can be written in a still more general form since 
there is no necessity for the interval of integration to be symmetrical 
about the origin. The only requirement is that the integral shall be 
taken over a complete period. In general the lower limit can have any 
arbitrary value t Q and the upper limit will then be / 0 + T. However in 
the majority of the cases which we shall be considering we shall use 
the form quoted above employing the symmetrical interval. 


1.5 Even and Odd Functions 

If f(t) is an even or an odd function then the series expansion has a 
special form. If fit) is even then by definition /(—/)= fit) and the 
expression fit)sinnt is odd. Hence the integrand for b n is odd. Now 
if we integrate an odd function over an interval which is symmetrical 
about the origin the result is zero, for if <f>it) is odd then 

jlj(‘)dt = j°_ dt+j T o 4,(1)di = J° <£(-/)<*(-/)+ \l<f>{i)dt 

= f r o <K-t)dt+j T o <Kt)dt 
and since <f>it) is odd <f>i-1) = — <£(/) 

/. jlj(t)dt = - jj <f>(t)dt + j T o <f>0)dt = 0 

Hence if the function fit) is even b n is zero for all n and we have 
simply a cosine series. Similarly if fit) is odd then a n is zero and we 
get a sine series. 
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Any function fit) can be expressed in terms of an even component 
fft) and an odd component f 0 (t) where 

m = £/*(')+/(-/)] 

and f 0 it) = WiO-X-t)] 

Thus in general we get a series containing both sines and cosines 
as already seen (section 1.2). 


1.6 Complex Form of the Series 

The occurrence of the sum of terms of the form cos«/ and sin nt 
suggests the possibility of deriving a complex form of the series using 
the relation 

cos 0+j sin 6 = exp ij6) 

We have 

oo 

f{t) = Ja 0 + S (a„cosnt+b n sinnt) 

H =1 


Consider the two terms a n cosnt+b„smnt. We can write 

a„ cos/!/+£„ sin nt = y/(p 2 n + b l) ^ 005 nt 

b„ . \ 

+- 77-5 — 

V(^+* n 2 ) ) 


V(p 2 n +bl)cos(nt-<f> n ) 


tan <t>n = - 

a* 


where 
Hence 

a*cos nt + b„ sin/if = W(p 2 n +*«)[exp {j(nt -<(>„)}+exp { -j(nt - <f> n )}] 

So far we have considered n to have only positive integral values. 
If we allow n formally to have negative values also then we have from 
the expressions in the last section 

a -n = &n* = -b n , and <t>_ n = -<f) n 
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so that 

a n cosnt + b n smnt = W(^l + ^)[expy{/;/ - <f> n } 
+exp y{(-«)/-A- B) }] 

Thus we can now write 


oo 

fit) = i^o+ S (a n cos nt+ b n sin///) 

n= 1 
00 

= i«o+ 2 W(al + b*)[expj{nt-<f>„} + expy{(-n)/-_ n) }] 

n= 1 
00 

= §a 0 + 2 4V(a»+ bl)expj{nt—(j> n } 

rt=l 

+ 2 4+*«) expy{( -n)l-<j> ( _ n) } 
n— 1 

= i fl o+ 2 W( a l+bl)expj{nt-<}>„} 

n— i 

+ 2 W(<*l+l>l)cxpj{nt-<f> n } 

n= — l 

= ^ + („?, + „ Ji,) ^ + eX PX W/ “ tn) } 

and noting that for n = 0 

i V +^2) expy(///— <f> n ) = since 6 0 = </> Q = 0 

we can write 

/W = 2 Witi+bl) expj(nl 


= 2 h/(ti + bl) exp ( -j<f>,) exp (jnt) 


00 

fO) = 2 C„ exp (jnt) 


or 
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where 


C n = W(a 2 „ + bl)zxp(-j<j>„) 

C n = i(a„—jb„) 


-IP 

2tt J_, 

-IP 


/(/){cos nt - j sin /»/} </f 


/(<) exp (-jnt) dt 

17 

Thus we finally obtain the relations 


/(O = 2 C„exp(y«/) 

n= —oo 


, if* 

"• 27rJ_, 


/(/)exp(-yVi/)A 


In this section we have allowed n to assume negative values in a 
purely formal fashion. A negative value of n in effect represents a 
negative frequency and we shall return later to the significance of this 
concept. 

As before we can extend this to the case of period T giving 
/(/) = jr C n expf/y«/j 

n— —oo ' 


i r‘‘+ T / 2tt \ 

c n = f(f) exp y-J—nt^dt 


where 
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1.7 Convergence of the Series 

In order to justify the Fourier series expansion we shall show that the 
infinite series 

00 

Ja 0 + 2 cos/*/+£„ sin wf 

n= 1 

does converge to the value/(/). Consider the finite series 


i^o+ a i cos t+b { sin 14-^2 cos 2 /+£ 2 sin 2 /+.. . + fl /J cosrt/+ 6 /l sin/f/ 
Define 


and 



/(/) cos ntdt 


b n 



f(t) sin ntdt 


where /(0 is defined in the interval -n^t^n and in this interval 
satisfies Dirichlet’s conditions. Outside this interval fit) is defined by 
the relation/(/+27 t) =/(f) so that the function is periodic and has 
period 2 xr. 

In the expressions defining a n and b m t is a dummy variable. Hence 
we can write 


a n cosnt+b n smnt 


J f(x) cos nx dx J cos nt 
+ |- J fix) sin nx dx j sin nt 


where/( jc) is now defined for - n < x < n and fix + 2 tt) = fix). This 
expression can be written 


ir 

"J-i 


fix){ cos nx cos nt +sin nx sin nt} dx 


=-r. 


fix) cos n(x—t)dx 


If S n it) denotes the sum of the finite series then 


S n it) 


.ir 


/(*){i+cos(*-/) + cos 2 (x-/)+... +cos nix-t)}dx 
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Let x-t = z 
and let 

T n = J + cos z+cos 2z+... + cos /iz 
2sin Jz.T n = sin $z+2 sin iz.cosz +... + 2sin £zcos nz 

= sin Jz—sin iz+sin f z—... —sin (/i — i) z+sin in + i)z 
= sin(rt+J)z 
. T = sin(/i+j)z 
n 2 sin iz 


Hence 


SniO 




„ N sin(/i + i)(*-/) 

fix) ^ -— dx 

2sini(x-/) 


We now investigate the behaviour of S n as n-+co. Consider the 
integral 

J%(£)sinoc idt, 

Suppose <£(£) is finite and satisfies Dirichlet’s conditions in the 
interval a ^t^b. When a becomes very large sina£ will go through 
a complete cycle with only a very small change in £ and hence in 
<£(£). If a becomes indefinitely large then <f>iQ can be considered 
constant whilst sina£ goes through one complete cycle and the 
integral will be zero, that is 


lim f 6 <£(£) sin a£ d£ = 0 

a—x© 


Returning to the expression for S n (t) and putting i(x-f) = C we 
have 

sin(/! + i)(x-/) 


l pi 

7rJ« 


2 sini(.x-/) 


-dx 


1 f ‘"" > 

Kw+O 2 suit 


Comparing this with the expression above ^^ corresponds to 
^(£). The only contribution to the integral (for large n) will be where 
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the condition <£(£) finite is violated, i.e. from the neighbourhood of 
points where sin £ = 0. The range of integration is from - i(^+ /) to 
\(7r—t) and is thus dependent upon the value of the parameter /. 
We wish to consider values of t in the interval - 77 < / < 77. Therefore 
the total range over which £ will vary is also - 77 to +tt and this 
includes three points where sin£ = 0, viz. £ = - tt, 0 , + tt. However 
we can exclude two of these points by considering the range - 7 t < t< tt 
which does not include the two points t = + 77 . Hence for this range 
there is only one contribution to the integral namely that from the 
neighbourhood of the point £ = 0. We shall deal with this case first 
and then complete the proof by considering the two points at the ends 
of the range of /, / = ± 77 separately. 

Since for the case - 77 < / < 77 the only contribution to the integral 
is from the neighbourhood of the point £ = 0 we can take the range 
of integration as - 8 to + 8 so that 


m- if 

TT J-a sm£ 

Furthermore in this interval we can put sin £ = £ 

s n (o = - r° At+2o- ia(2n+iH dc 

* J-B £ 


-if 

” Jo 




Now as £->0 from negative values /(/+2£) becomes /(/_), the 
limit of the function f(t) at the point t when approached from the 
negative side. Putting (2n+ 1 )£ = f and letting £-*0 


if 


A/+ 20 - (2 ^ +1)C rfC 

* £ 


- ;A>-) f ^ 

n J-(2n+1)a C 


It 


As n -> 00 this becomes 
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Similarly it may be shown that 


if 

^ Jo 


'/»+*> 


Hence we have finally 

lim S n (t) = i{/0+)+/0-)} 


If the function /(/) is continuous then /(/+) =/(/_) and S n {t) ->/(/). 
At points of discontinuity S n tends to the average value. 

The expansion has been justified for — 77 < t < tt. To complete the 
proof we must investigate the points t = ± 77. We require the values 
lim S n ( ± 77). Consider first the value for t = + 77. We have 

n—*co 




sin (n + j)(x-t) 
2 sin i(x—t) 


Substituting i(x—t) = — £ 


1 f 

*n(‘) = - 

77 Ji 


1 


W + v) 
1 /•!('+*) 


At-20 


1 f* ( '+”> sii 

At- 2D- 

"JUt-rr) 

_, . If” ^si 
S„(v) = - /(*r-2£)- . r 

■n Jo sm £ 


sin ( 2 /i+l)£ 
sin £ 

sin( 2 /«-t-l)£ 


sin £ 
sin (2 n+ 1 ) £ 


</£ 


(~d0 

di 


Following the same procedure as before the only contribution to 
the integral will be where sin £ = 0 , i.e. in the neighbourhood of 
£ = 0 and £ = it. 


A SJfr) = - f 

* J 0 


V(.-20S!M« 

sin£ 


1 r 
w j»~ 


y sm£ 






















22 


AN INTRODUCTION TO FOURIER ANALYSIS 


Substituting £' = tt+ £ in the second integral 


SJM 


- - f /("“ 

^ Jo 


2 D 


sin ( 2 /i+l)£ 


sin£ 




-- f /(- 7 T- 2 D: 

j-* 


, sin( 2 /i+ 1 )£' ,, 

> *r * 




'dt, 


As £ —>0 through positive values /(w-2£) becomes/(w_) and 
sin £->■£. As £'-*•() through negative values /(-«•—2{') becomes 
/(-w + ) and sin 

1 r , N f a sin( 2 «+l)£ 

■ ; /<0 J. —r 

,I„ , f sin( 2 »+ 

+ ; /( ~’ J J-r —r—* 

Finally putting (2 a+ 1)£ = f and (2 /j +1)£' = £' and letting /i->oo 
we obtain 

S„(n) = ^/(ir_) \f{ “ w + ) J ^ ' 


= *{/(»_)+/(-»+)} 

But since/(t+2ir) =/(/), /(- 2 r + ) =/(tt + ) 

£«(*) = i{/( 7r +)+/( 7r -)} 

Similarly it can be shown that 

SJL-") = i{/( -Tr +)+/( _7r -)} 


1.8 Multiplication Theorem 

Using the complex form of the series we have 

fit) = 2 C„exp(jnt) 

n= —oo 

C n = ^ J /(/) exp (- 7 / 1 /) dt 
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Consider the average value of the product/, (t) ./ 2 (/) over a complete 
period. This is given by 

^ J_ fl(t)f 2 U)dt = ^ J | ^ ( C l)n eX PC/«0|/ 2 (0* 

= 2 ( C, )«{^ J_ / 2 W exp (/«/)<*| 

= £ (C,UC2)_„ 

n— —00 

i.e. 

^ r /,(o/ 2 « ^ = 2 ( c «)»( c 2 )-» 

" v n= -00 

Thus the average value of the product of two functions over a 
complete period is equal to the sum of the products of the coefficients 
(Ci)n(Q) -n where the C’s are the Fourier coefficients in the complex 
series representation of the functions. 


1.9 Parseval’s Theorem 

From the multiplication theorem deduced in the previous section we 
have by putting/, =/ 2 =/ 

00 

i jm’*- 2 c * c - 

** w #i= —00 

Now from section 1.6 we have 

C„ = K«n-A) 

Also we have a_ n = a n and = —b n 
Hence 

C = C* 

where C* denotes the complex conjugate. 

c n c_„ = c„c* = |c„| 2 
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2tt 



CO 


2 i c -i’ 

1 = —00 


Writing C n in terms of the components a m b n 

\ c „\ 2 = Itf+bl) 

••• S |Q| 2 = s 2 i(«M) 

fl=-oo n=—co n= 1 

hr f" (W* = i fl 0 + 2 

A7T J-V „= 1 

This is Parseval’s theorem. 


1.10 Power in Complex Waveform 

Consider the term a n cosnt. The power in this waveform is given by 
l C n 

^ I ( a„cosnt) 2 dt 

For example if a n represents a voltage then the expression gives the 
power which would be developed in a resistance of unit value. 

Now 

l C v 

-jJa n cos nt y dl = ia 2 n 

Similarly the power in the component 6 „sin//f is given by \b 2 n . The 
power in the complex waveform /(/) is given by 

hfiwr* 

which by Parseval’s theorem is equal to 

ial+ 2 Itf+b*) 

n— 1 

Hence the power in the complex waveform is equal to the sum of 
the powers in the Fourier components. 


CHAPTER II 


Analysis of Periodic Waveforms 


2.1 Introduction 

We have seen in Chapter I that a function/(/) which is periodic, i.e. 
which is such that f(t+T)=f(t) can be represented in terms of a 
Fourier series 


/(/) = 



cos 


2tt , 2tt \ 

— nt+b n sm — nt I 
T T ) 


In this chapter we shall apply these ideas to analyse some of the 
waveforms commonly occurring in the physical sciences. 


2.2 The Square Wave 

A sketch of this waveform is shown in Fig. 2.1. It is assumed that the 
amplitude of the wave is unity and the period T. With the choice of 
origin shown in Fig. 2.1 the function is anti-symmetric or odd. 




+ 1 






- 


0 1 
2 

T 



-l 



Fig. 2.1 

The square wave—anti-sym¬ 
metric form 


Hence the Fourier expansion will be a sine series and we have a n = 0 
for all n (including n = 0). Using the expressions developed in Chapter 
I we have 
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Furthermore /(/) is defined by the expressions 

m = -i, -\t< t < o 

= 1, 0 < t< iT 

and f(t+T)=f(t) 


Therefore 


2 f f 0 2rr f r/2 h \ 

= p|J ^(-Osin—;»/<*+I (+ I)sin— 

T/2 


-r /2 

2/r r 2 W 
= rfc cos 7'" 


1° IT 2tt f 

J-r/2 l>* r Jo 


= -{[1 — COS(— //7r)]+[l — COS(/27r)l} 

7I7T 


= -(1—COS/27T) 

JllT 


.*. b n = 0, n even 
= —, n odd 

nir 



Fig. 2.2 

The square wave—symmetric form 


We could equally well have chosen the origin as shown in Fig. 2.2. 
In this case the function fit) would be defined as follows 

/(f) = -1, -i T<t<-\T 

= +1, -±T< t < i T 

=-l, jrc/cjr 

and /(/+D = /(/) 
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The function is now even and hence b n = 0. We have 
rr /2 


27 


«o 


mdt 


2 1 r ~ T/4 r w r 7 */ 2 

= J (-n<*+ (+d</h- (-n* 

1 U-r /2 J-r/4 Jr/4 


= 0 

It can be seen that this must be the case since ±a 0 represents the 
average value of fit) over one complete cycle (the d.c. term in electrical 
terminology) which is clearly zero. The rest of the coefficients are 
given by 

r*V2 2 t r 

/(/) cos —nt dt 
-Til T 


2 rn 
° n Tj-j 

2 r-r/4 2ir 2 f r/4 27 t 

= t\ (-l)cos—»/<//+- I (+l)cos—«/<// 

1 J-T /2 1 1 J-r/4 1 

2 f r ' 2 , ts 2n J 
+ ~ (-l)cos —ntdt 
1 Jr/4 * 


= -(- 

T. 2 n } 

,-r/4 j 

r. 2 n 7' 4 

r. 2 n 

n 


- sin—«/ 

+ 


sin —nt 

i 

7 rn\ 

L T \ 

1—r/2 1 

L T J-r/4 

L T . 

lr/4/ 


If . / 7i7r\ . //J7r\ . / mt\ . /nn\\ 

-™r s,n (-T) +s " , (y)- s "’(-Tj +s,n (Tj) 


4 . 727T 
= —sin— 
irn 2 


= 0, /i even 


= + - 


Mr 


n = 1,5,9,... or 7i = 4&-3 


=-, 71 = 3,7,11,... or n — 4k—\ 

mr 


where k= 1,2,3,... 
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By a suitable choice of origin we can expand the function either as 
a cosine series or as a sine series. The origin can of course be chosen 
anywhere yielding in general a series containing both sines and 
cosines. 

2.3 Waveform Symmetry 

In the previous section we deduced that for the square wave the even 
harmonics are all zero. It is interesting to investigate the condition 



Fig. 2.3 

The effects of Harmonic content (a) Odd harmonics; (b) even 
harmonics 


for which this property occurs a little further. We can readily sketch 
the effect of the second harmonic and the third harmonic on the 
waveform of the fundamental. This is shown in Fig. 2.3. These two 
waveforms have a different type of symmetry. For the third harmonic 
case (Fig. 2.3b) the reflection of the second half of the waveform, (that 
between \T and T), in the t axis produces a waveform (shown dotted) 
of exactly the same shape as that between 0 and £ T. In the case of the 
second harmonic however the two corresponding waveforms are 
mirror images. We can express the type of symmetry exhibited in 
Fig. 2.3b by the relation 


fo+m = -/(o 


^1 K> 
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We then have 
2 f 772 


2 r T 

a "" rj_ 


f(t) cos ^ntdt 
r/2 T 

2 n 2 /*r/2 2 

/(/) cos — ntdt+— fit) cos-^ ntdt 

r /2 T T Jo T 


2 C ri2 2n 2 r*V 2 

- t Jo ^ t+ ^ cos Y n ^ + ^ dt+ fj 0 

2 r m /2ir \ 2 r T < 2 2IT 

= ^J 0 -/(Ocos^— nt+nirjdt+- I fit) cos —ntdt 

= ^jo ^( cos T'“ ~ °° s (^f nt +W,T |J dt 

2 C T ' 2 2n 

= j-I /(/) cos — w/(l — cos rnr) dt 

= 0, n even 

Similarly 

2 C m 2n 

'« = fit) sin —ntdt 

1 J— 772 T 

2 r° 2n 2 r r/2 2n 

= - fit) sin — ntdt + - fit) sin - ntdt 

1 J—r/2 T T Jo T 

2 f r ' 2 s 2ir , 2 C TI2 2n 

~ f Jo /(/ + £T) sin —nit + iT) rf/ + - I /(/)sin — ntdt 

2 f r ' 2 , /v . (2-n \ 2 C TI2 2i t 

~ TJo f(0sinl—nt+nnjdt+- I fit) sin—nt dt 

2 rv2 | 2 ir . / 2 it \) 

= ^.J /(Ojsin — nt-s\n\—nt+nir\\dt 

i 


r /2 2tt 

f(t) sin — nt (1 - cos mr) dt 


Jo T 

= 0 , n even 
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Thus we see that waves exhibiting symmetry expressed by the 
relation/(/+ \T) = —fit) have no even harmonic content. The push- 
pull amplifier is an example of a practical circuit which produces 
waveforms with no even harmonic distortion. 

2.4 The Triangular Wave 

A sketch of this waveform is shown in Fig. 2.4. This can be expressed 
analytically as 

fit) = i+p, -ir</<o 

4 

= 0<t<iT 



Fio. 2.4 

The triangular wave 


We see that/(- /) =/(/), i.e. the function is even and hence b n = 0. 

- ■ - KJ'+Phiri'-P)* - ° 


This too could have been predicted since Ja 0 represents the average 
value of fit) over one complete cycle (the d.c. term in electrical 
terminology) which is clearly zero. 




2 f 772 2«r 

- fit) cos ~ nt dt 

T J-TI2 T 

2 r T l 2 2ir 2 |*° 4 2ir 

?J- ra cos 7" , ‘" + rJ. TO f-'“ s r“‘* 


2 C T ' 2 4 2n J 
-- — / cos —nt dt 

TJ o T T 
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8 |*° 2ir . . 8 r T ' 2 2 it 

°n = - 2 j + r2 (-/) cos -«^(-/)-- 2 | o t COS—It tdt 


16 

C Tl2 2tt 

T 2 , 

tcos-—ntdt 
Jo T 

161 

T T . 2tt 1 r ' 2 

= — — < 

\t -sin— nt\ 

r 2 ( 

1 2 t r/i T Jo 

8 r 

*3 

_t 

5 

KJ 

■ z=— • 

- -— cos — nt 

Tirni 

2nn T Jo 


= ^(l-COS/ITT) 


n 2 n 2 ' 


n odd 


= 0, n even 

Note that this function satisfies the relation /(/+ JT) = -fit) and 
hence all the even harmonics are zero as we showed generally in the 
last section. 


2.5 The Saw-Tooth Wave 

The saw-tooth is a waveform which occurs quite frequently in elec¬ 
tronic engineering. It is sketched in Fig. 2.5, and can be represented 
by the function fit) defined by 

fit) = p -IT< 0 < IT 



Fio. 2.5 

The Sawtooth wave 
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In this case we have an odd function since /(—/) = -/(/). Hence 
we have a n = 0 and b n given by 




2 f 772 2 t r 4 r 


r » . 2 » , 
/sin —ntdt 
r/2 7 1 


4 (I 

f T 2ir 

1772 T 

/*772 

2tt \ 

~A\ 

— - — tcos—nt 
[Inn T 

l-r /2 2irn „ 

l-r /2 

cos — ntdt\ 
T 1 


2 fr 


7777*1 2 


- cos (/?7r) 4- - COS (— W7r) 


=-COS (//7r) 


mr 


2 

.. — + , 
mr 

n odd 

2 

tnr 

n even 


2.6 The Infinite Train of Pulses 

This is another commonly occurring waveform in practical tech¬ 
nology. We shall first assume pulses of unit height, duration 2 r and 
period T. 

The expression for/(/) is 

/(/) = 0 -i r</<-r 

= 1 —T < / < T 


= 0 T < t < \T 

/(/+r)=/(/) 


The function/(/) is even and hence = 0. For the case shown in 
Fig. 2.6 there is a d.c. term and this will be given by \a Q where 



cm 

J-r/2 


Mdt = | 
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The rest of the coefficients a n are given by 


2 f 772 „ x 2 tt . 2 f 7 2 tt 

= f]_ Tj /( t)cos Y ntd * = t cos Y ntdt 
1 r. 2tt r 
mrl T J_ t 


2 . 27t 

.*. a„ = — sin— ht 
mr T 


It is interesting to consider the limiting case of very narrow pulses, 
or impulses. Suppose we consider pulses of unit height and duration 


n n n 




-ror 


H&- 


Fig. 2.6 
A train of pulses 

2 r and allow r to become very small. In this case we can write 

2 . 2t r 2 27r 4r 

o = —sin — m — — — nr = — 
n mr T mr T T 

and the magnitude of the a n become independent of n. Thus the 
spectrum consists of an infinite series of lines of constant height and 
spacing. The original function of time we considered had exactly the 
same form. Functions of this type which have the same form when 
they are represented as a function of time as they have when repre¬ 
sented as a function of frequency (the spectrum) are said to be 
self-reciprocal. 

It may be noted that when r->0, a n ->0 and all the spectra have 
zero amplitude. This is because the original time waveform which is 
composed of spikes of finite amplitude is now represented as the sum 
of an infinite number of spectral terms all of which have become 
equally important. We can also look at this in terms of energy. If the 

B 
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original time waveform consists of spikes of finite amplitude and zero 
width then the energy in the waveform must be zero. We know 
furthermore from section 1.10 that this must be equal to the sum of 
the energies in the spectral components and hence these must all be 
zero. 

We could form a train of unit impulses by defining the amplitude of 
the pulses as 1/2r instead of unity. The area under each pulse would 
then always be one. The spectral components will then have value 

1 2 . 2 tt 

a n = --sin — nr 

2r IlTT T 

which as t-*0 becomes 

^ 4t _ 2 
a "~ 2r T ~ T 

The infinite series of impulses can therefore be written 
x 1 2 2tt 

m = f + fZ cos f nt 
"- 1 

Under these conditions the amplitude of the original pulses becomes 
infinite, the energy per pulse remains finite, and since there are an 
infinite number of pulses in the train the total energy in the waveform 
is infinite also. In the frequency representation the amplitude of the 
individual components is finite, representing finite energy but again 
there is an infinite number of them so that the total energy is infinite. 

2.7 The Rectifier 

Rectifiers are used in electrical engineering to convert alternating 
currents and voltages into direct currents and voltages. The rectifier 
conducts only when the voltage is applied in one direction and not 
when the polarity is reversed. Thus if the applied voltage is sinusoidal 
the rectifier conducts only on the positive half cycles. It is common 
practice to use a combination of two or more rectifiers to produce 
what is known as full-wave rectification. In this scheme the negative 
half cycles are effectively reversed to produce a waveform as shown in 
Fig. 2.7. 
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Note that the period of the cosine wave from which the waveform 
of Fig. 2.7 is produced is 2T Thus it is seen that the fundamental 
component of the rectified cosine wave is double that of the original 
cosine wave. (In electrical terminology the ‘hum’ produced by full- 
wave rectifying the 50 c/s supply has a frequency of 100 c/s). 



-£r o l t 
Fig. 2.7 

A rectified cosine wave 
The function representing this waveform is 

/(/) = cos“ -iT< t < IT 

The d.c. component is given by \a Q where 


2 r* it 2T . tt 7/■ 

““ ■ rJ- ro “V‘" ‘ 


1 772 
Til 


4 

a Q — ~ 

7T 

Thus the d.c. term is I/tt. 

The function is even and therefore b n = 0. The harmonics are given 
by a n where 


1 r^ii f ir 

= - J |cos(2u—l)-/+cos(2/i 


?)■" 


| r /2 
772 
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If n is even 




and if n is odd 


2r_i_i_i = _ 

7 q_ 2 /i — 1 2/i+ 1J + 7r(2n — 


rr{2n— l)(2/i + I) 
4 


D( 2 «+l) 


2.8 Gibbs’ Phenomenon 

The series developed above represents the function /(/) only if an 
infinite number of terms are taken. We shall investigate now the form 
of the curve when only a finite number of terms are taken. We can 
illustrate the phenomenon discovered by Gibbs by considering the 
square wave we investigated earlier (section 2 . 2 ). 

For simplicity consider the square wave of unit amplitude and 
period 2n f i.e. 

/(/) = -1 -7r < / < 0 

= +1 0 < / < 7T 


Consider the sum of a finite number of terms of the series. From 
section 1 .7 this is given by 






sin h(x-t) 

Hence for the square wave we are considering 


1 f w si: 

w “si- 


sinQt+Dfrf-Q , j fo sin («+$)(*'-/) 

sin i(x' - /) 


'-hL 


l dx* 


jo sinJOt—0 
Substituting * — / = 0 and x' — t = — 0 r we have 

i„s M - r*siss»** + r *<(,,+»* 

J-t J w+ , sin 

Both 6 and O' are dummy variables in this expression. Using 0 as 
the dummy variable throughout and rearranging we have 


2 wSM = P S ± (n+ J‘ )e de+ r 

J-, sini^ J„ +l sin}0 


dO 
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Now we know that when n becomes very large S n (t) except at 
the points of discontinuity, i.e. in this case at / = 0 , when 

S n (t) ^ i{/(°+)+/(°-)} 

We shall now investigate how S„(/) behaves in the neighbourhood of 
the discontinuity, i.e. t = 0 , when n is finite. 

The first integral, for t small, is evaluated over the region where 
0 = 0. At this point both the numerator and the denominator of the 
integrand become zero. The value of the integrand at 0 = 0 is obtained 
by taking the ratio of the derivatives of numerator and denominator, 
i.e. 

[{(« + i) cos ( w + i) 0 }/i cos = 2/1+1 

The second integral is evaluated over the region where 0 = w. Here 
the magnitude of the integrand is approximately unity. Thus for all 
but the smallest values of n we can neglect the contribution from the 
second integral in comparison with that from the first. Thus for the 
conditions we are interested in investigating we can write 


2irS„(t) 


f' sin(j 

~ J-» si 


(»+*)« 


sini# 

Since the integrand is even in 0 we have 

C 1 sin (n +1)9 

2irS n {t) = 2 

Jo sinjd 


d0 


C' sin(i 

2wS » ( ' ) = 2 Jo^ 


dO 


and since we are interested in the region near the origin we can 
replace sin $0 by id giving 

‘ f sin(/i +£)0 


2t tS, 

Substituting n + J = m 
7 tS, 


r si 

■.0) = 2p 


ie 

sin m0 


-d0 


d0 


This is related to the sine integral, an important integral defined by 
the expression 
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This integral is obviously zero when f = 0 and when f->-oo can be 
evaluated by the method of contour integration yielding a value in. 
At other points it can be evaluated by numerical, methods and has 
been tabulated. The form of the curve is shown in Fig. 2.8. 
Considering now the form we deduced above 

,. f' sin mO . 

-i~r M 

Substitute mO = k 

C mt sink lt 

g(t) = I —£-dk = Si(mt) 



Fio. 2.8. 

The Function Si(t) 


Thus as m (or n ) increases the curve is compressed into the origin 
but the general form of the curve does not change. 

Returning now to the function we were considering, the finite series 
approximation to it is given by 


Thus 

and as /j — ► oo 


S„0) = 


2 sin A: 

Jo k 


S„( 0 ) = 
S„0) = 




1 


Hence at the origin the value of S n is zero. It then rises rapidly as t 
increases and overshoots the value unity, the value of /(f) in this 
region, passes through a maximum and oscillates about the line 
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/(/) = 1 with decreasing amplitude. The value of the overshoot is not 
dependent on n. As n increases the position at which the maximum 
occurs moves nearer to the point of discontinuity. If we consider a 
point removed from the discontinuity then as we have already 
mentioned as n increases the curve of Fig. 2.8 is compressed into the 
origin and therefore at the point under consideration the smaller the 
amplitude of the oscillation and the closer the value of S n approaches 
one, the value of/(/). 

We can see more generally that S n will produce an overshoot in the 
neighbourhood of a discontinuity as follows. Consider a function/(f) 
which has a discontinuity at f = 0. Suppose for simplicity the function 
is odd. Let/(0 + ) = A. We then have/(0_) = - A and/(0) = 0. Let 
/(f) further be periodic with period In. The function /(f) can be 
expressed as a sine series 

00 

fit) = 2 b r s\nrt 

r-1 

Hence 

/(f)—5 n (f) = £ b r sin rt 

r««+l 

Now 

00 - 

f” {/(f)-S„(f)}sinfrff = S b r sin rt sin tdt=0 

Since/(f), S„(f) and sin f are all odd the integrand on the left hand 
side is even and hence 

f {/(f)— S n (l)} sin f = 0 

Over the range of integration in this expression sinf is positive. 
Now the function S„(f) is the sum of a finite number of sinusoidal 
terms and hence 5 n (r) will approach zero when / is small. The function 
fit) on the other hand jumps discontinuously to a finite value. Hence 
for small positive values of /, fit) - S(t) < 0 and will therefore give a 
negative contribution to the integral. This must be balanced by a 
positive contribution since the value of the integral is zero. Hence 
S n (t) must overshoot the value of fit) in some other region of the 
interval. 












CHAPTER III 


Fourier Integrals 


3.1 Introduction 

We saw in Chapter I that a function fit) defined over the region 
— i t ^ IT could be represented as an infinite series 


i°o + 2 {a n cos y nt+b„ sin y n/j 


or letting 2ir/T= to 


CO 

ia Q + 2 cossin//<*>/) 

H"1 


We can plot a„ and b n as functions of frequency to. These plots 
consist of a number of discrete spectral lines occurring at a constant 
spacing of to = 2n/T. As the region over which the function is defined 
becomes larger so the spacing of the lines becomes smaller. In the 
limit when the function is defined for all t, T-vooand the spacing of 
the spectral lines becomes zero. Thus the plots of a n and b„ are no 
longer discrete but become continuous and the sum becomes an 
integral. We shall now use this idea to derive the Fourier Integral 
Theorem as a limiting case of the series. 


3.2 The Fourier Integral Theorem 

If we consider a function/(/) defined in the region -n/Q^tsS. n/Q 
and obeying Dirichlet’s conditions then we have already seen that we 
can represent/(/) by the series 

00 

f(f) = S (Pn cos nQt + sin nQt) 



J ir/a 

fit) dt 

-nfS3 

40 


where 
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q 

"■'ri- 


Hence 


Q C n,a 

n/O 


fit) cos nQt dt 
fit) sin nQt dt 


O Mn Ot f"l Q 

fit) = — /(A) / — Jcos nQt /(A) cos nQX dX 

4- sin nQt f /(A) sin /z^A Ja| 

J-ir/0 I 

h . 2 ri- 


_ Q 

~ 2tt J_ 


fiX)dX+ 


/(A) cos nQ(t-\)d\ 


where A is a dummy variable. 

If we now take the limiting case f2-*0 this expression becomes 
valid for all t . The spacing of the spectral lines in the plot against co is 
Q and hence as Q ->0 the lines become more and more closely spaced. 
In the limit the spacing can be represented by dco the increment of 
the continuous variable a>, and we get a continuous spectrum. Now 
n assumes all integral values from one to infinity and hence as Q ->0, 
nQ varies continuously from 0 to oo (provided we make and 
n -> oo in such a way that nQ -*oo). We can therefore represent nQ by 
the continuous variable w. The summation in the limit passes over to 
an integral and the limits as we have just seen will be 0 tooo. Hence the 
second term in the expression can be written 



/(A) cos <oit-X)d\ 


The integral in the first term becomes /(A) dA in the limit and 
this is finite since /(A) satisfies Dirichlet’s conditions. Hence as Q-+0 
the first term in the expression for/(/) also tends to zero. Thus we 
have 

/(/) = - I dco I /(A) cos co(t— A) dX 

TT Jo J—00 

This is Fourier’s Integral Theorem. 

B* 
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3.3 Proof of the Integral Theorem 

We shall now prove the Integral Theorem by a procedure analogous 
to that used in investigating the convergence of Fourier series. 
Consider the integral 

Ia(!) = JJ* dto /(A) cos <u(t - A) dX 

Reversing the order of integration we have 

IaO) = f ^ / 0 ) ^ J Q cos <i)(t— A) dco 

- 

. r /( a>!“I^->,» 

J — co t—A 

We now let A-> oo. Using again the theorem 
lira |%(£)sina£d£ = 0 

a-*oo 

provided </>(£) is finite in the interval ab we see that the only contri¬ 
bution to the integral will be in the neighbourhood of the point where 
/ = A. Substituting A - / = 6 we have 




f(f+0)—^-d0 


As /4->oo the only contribution to this integral will be from the 
region near 0 = 0 . 


>*-►00 J-0 V 


sin/40 


-£ 


/ <,+9 >. 


1’ 


sin/40 


d 0 + | /(r+ 0 ) 

sin/40 


d6 


d6 
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Writing A6 = (f> we have 

A —► oo J-oo 9 Jo 9 

. . f° sini . . f® sin<£ ,, 

= /(*-) J oo "^-^+/( / +)J 0 


= 5 /(' + )+/('->} 


Hence 


| /*oo /*oo 

- \ da>\ 

»J 0 J—O 


/(A) cos <o(/—A) dX = *{/('+)+/('-)} 


If f(j+) = /(/_), i.c./(/) is continuous at this point the expression 
becomes equal to f(t) so that 


i /•oo /*oo 

71 JO J-a 


f(A)cosw(t-A)dA = /(/) 


3.4 Sine and Cosine Transforms 
We have from the integral theorem 

| /’oo /’oo 

/« = - 

77 J 0 J-« 

Expanding the cosine 

30 

/(A){cos a>/cos cuX + sin wt sin to A} dX 


/(A) cos w(t-X)dX 


1 /’oo /’oo 

/(/) = - t/oA 

77 Jo J-« 

i p* 

.%/» = - 

77 Jo 


Jo 

+- 


/•oo 


/(A) cos cuA JA 


cos cot dco 

] /’oo /•oo 

-- I sin cot dco I /(A)sincoA</A 
77 JO J-oo 

We define the Cosine Transform a(co) as 

a(o>) = /(A) cos o>A dA 
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and the Sine Transform b(<o) as 

b(<o) = j ^/(A)sin coXdX 

The function/(/) can now be written in terms of these transforms. 


1 f 00 

/( t) = - {a(aj)cosa>/+^(cu)sina;/}^cj 
" Jo 

The integrand contains the sum of a term in cos cot and sinw/ so 
this can also be written in the amplitude-phase form 


where 

and 


1 f°° 

fit) = - I R(<o) cos [u)t+</)((*))} d<D 

77 Jo 


Rico) = V[{a(co)} 2 +{b(<o)} 2 ] 


tan (f>(co) = — 


bjco) 

a(w) 


Note that the region over which we integrate a(co) and b(<o) is from 
0 to ooand hence we are only concerned so far with the case <o > 0 , i.e. 
positive frequencies. In the next section we shall need to extend the 
range of (o to the negative values. The significance of negative 
frequencies will be discussed later. 


3.5 Complex Transform 

Using the results of the previous section we have 

fit) = - J \/[{«M} 2 +{^M} 2 ]cos {<ot+</)(<o)}d<o 

= ^ J o V'[{ a (“»)} 2 +{*(«)} 2 ] i[exp{yaj/+y^(to)} 

+ exp { —jwt -j<f>(w)}] dw 

We now formally allow <o to assume negative values and define 
a(<o) and b(<o) as in the previous section for all <o. We note that a(co) 
is even and b(co) and </>(co) are odd, i.e. 

a(— (o) = a((o) 
b(—oj) = — b(oj) 

<f>(-<o) = <f*(to) 
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Hence we can write 

/(/) = ^ f V(a 2 + 6 2 )[exp (jwt) exp {j<f>(w)} 

+exp (-you/) exp { -j<f>(u )}] dw 

= ^- f \Z{a 2 + 2 > 2 )[exp (Jwt) exp {j<f>(w)} 

2ir Jo 

+ exp {y( - co)/} exp {j<f>( - a>)}] dw 
1 f°° 

= — \/(a 2 + b 2 ) exp (jwt) exp {j<f>(<o)} d<o 

2tt Jo 

+ T" f V(* 2 + b 1 ) exp (J<*>t) exp {j<f>(co)} d( - <o) 

Jo 

where (o has been replaced by — <o in the second integral. Hence 

i r°° 

/(f) = — V (a 1 + b 2 ) exp {jwt) exp {j<f>(oj)} dw 

2n Jo 

i r° 

+ — V( a2 +b 2 ) exp (Jwt) exp {j<f>(w)} dw 

2tt J —oo 

.*./(/) = ^- I V( a2 +b 2 )txp {j4>{w )} exp (Jwt) dw 
2tt J_oo 


Now define 


Then 


F(w) = \/[{a(a >)} 2 + {ftCto)} 2 ] exp {j<f>(w)} 


At) 


l_ P° 

“ 2n J_a 


F(co) exp (jcoi) dto 


The expression for F(co) can be simplified 

F(w) = -\/(a 2 +b 2 ) exp(y^) = V(° 2 + b 2 ){cos <£ +y sin 

= V ( a 2 + * 2 )[ A /( a 2 + -J ^ a 2 + 6 2 )] 

.*. F(oj) = a(co) —jb(co) 
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Fiw) can also be expressed in terms of/(/) 

F( w) = a(co) —jb(aj) 

= I fit) cos wt dt-j f°° fit) sin wt dt 

J— oo J —oo 

F(co) = /(/) exp (-yW) dt 

The function Fiw) is known as the complex Fourier transform of 
fit) and /(/) is the inverse transform of Fiw). The two relations are 
expressed by the formulae 




F(cd) expijwt)dw 


fit) exp ( -jwt) dt 

The two functions/(/) and F(a>) are known as a Fourier pair. 

The Fourier transforms are special cases of the general integral 
transform/(a;) of/(f) (see for example Tranter, 1951), defined by 

/(«) = 

The function K(w,t) is called the kernel of the transform. In the 
case of the cosine transform Kiw,t) = coswt, and the sine transform 
K(a>, 1) = sino>/. In both cases a — 0 and b = oo, so that 


/(<o) = J^° fit) cos wt dt 
f(aj) = fit)s\nwtdt 


Cosine Transform 
Sine Transform 


For the complex Fourier transform Ki<o , t) = exp(- ja>t ) 9 a = - oo 
and b = a> 


fi<x)) = J ^ /(/) exp ( —jcjt) dt Complex Transform 

In general in this book we shall be concerned with functions fit) 
which are real functions of time. In this case F(co) will in general be a 
complex function of the variable co (which is real) and can be written 
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in the form a(w) - jb( a>), where a and b are both real functions, as we 
have shown above. However there are cases of interest in which/(/) 
may be complex; for example in modulation theory where we use the 
function exp ijwt). hit), where h is real. 

Since F(a>) is in general complex, if we wish to plot F(o>) as a 
function of frequency we must plot two graphs, one representing the 
real part and one representing the imaginary part, i.e. we plot aioj) 
and biw). Note however that biw) is not the imaginary part of F(o>) 
since F(o>) = aia>) - jbiw). Thus - biw) is the imaginary part of F(w). 

We shall now consider some of the properties of F( w). If fit) is an 
even function then/(/)sinojf is odd in t. Hence 

biw) = j™^ fit) sin wtdt = 0 

Thus for the case of an even function the complex transform F(co) 
degenerates to a pure real function. Similarly if fit) is odd then 
aiw) = 0 and so for an odd function the transform F(o>) becomes pure 
imaginary. 

We have already noted that aiw) is even and biw) is odd. Hence if 
fit) is a real function so that aiw) and biw) are real then we have 

F(-o>) = ai-w)-jbi-w) = aiw)+jbiw) = F*(co) 

so that 

F*iw) = Fi-w) 

where F*iw) denotes the complex conjugate of Fiw). 

If fit) is complex then let 

fit) = hit)+jk«) 

where hit) and kit) are real. 

Fiw) = ^fit)cxvi-jwt)dt 

= J°°^ {/)(/) 4- jkit)}(cos wt - j sin wt) dt 
= [{/*(/) cos wt + kit) sin wt} 

+ j{kit) cos wt - hit) sin co/}] dt 
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F(-co) = J 00 ^ [{/*(/) cos co/—&(/)sinco/} 

+ j{K0 sin cot + k(t) cos co/}] dt 
F*(a>) = j ^ [{/*(/) cos cot+ k(t) sin co/} 

+ j{h(t) sin co/ - k(t) cos co/}] dt 
Hence if F(— co) = F*(co) then we have the two conditions 

I k(t) sin cot dt = — f°° k(t) sin cotdt 

J— 00 J—00 

J 00 ^ k(t) sin cotdt = 0 

/. K-t) = k{t) 


and 


I k(t) cos cotdt = - I k(t) cos cotdt 

J —00 J —00 

/. J°° k(t) cos cotdt = 0 

/. k(-t) = -*(/) 

Hence from these two conditions we have 


k(-t) = k(t) = —kit) 
k(t) = 0 

Thus if F(-co) = F*(<o) then k(t) = 0 and/(/) is a real function. 
We therefore have the result that F(-co) = F*(co) if and only if /(/) 
is a real function of /. 


3.6 The Concept of Negative Frequencies 

In the previous section we formally allowed co to take negative values 
in order to develop the formula for the complex transform. A similar 
procedure was also used in the derivation of the complex form of the 
Fourier series. We shall now examine the concept of negative fre¬ 
quency a little more closely. 
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Since in each case the need for negative frequencies arose when we 
developed the complex form of the formula we shall start by discussing 
the complex representation of a sinusoidal type of variation, viz. 
exp0’o>0- This is a complex function and can be written in terms of 
its components thus 

exp 0*^0 = coscj/+ysinco/ 

This is a well-known relation. We can plot this on an Argand 
diagram in which the real part of the function is plotted along the 
x axis and the imaginary part along the y axis. The form of the 
function expO’^0 shows that the amplitude is constant and its phase 
<£(/) is equal to cot, i.e. it increases linearly at a rate determined by co. 
Thus the function represents a point which moves round a circle in 
the Argand diagram at a constant rate in the positive direction of 
rotation, i.e. by convention counter-clockwise. It is now fairly easy 
to attach an interpretation to a negative frequency. Thus the phase 
of the function exp ( - jcot ) is <£(/) = - cot, i.e. it increases at a constant 
rate but in the negative direction of rotation (clockwise). We some¬ 
times think of a quantity which is specified by an amplitude and a 
phase as a vector. Hence we have the idea of representing a sinusoidal 
variation in terms of a rotating vector. If we require a real co¬ 
sinusoidal variation then we add together two vectors rotating at 
the same rate but in opposite directions, i.e. a positive and a negative 
frequency. This is represented analytically by the well-known relation 

cos cot = £{exp (Jcot) + exp (-jcot)} 

The complex spectrum of cos cot is pure real (we saw above that 
this is true for all even functions) and consists of two spikes one at 
+ co and one in the negative range of frequencies at — co. The spectrum 
is a line spectrum because coso>/ is periodic. It is easily seen that the 
magnitude of the spectral lines is tt since 

1 f 00 

fit) = — F(co) exp (Jcot) dco 

2“3T J —00 

= ^-[77 eXP (jo) Q t) + TT ex P ( —jcO 0 f)] 

27 T 

= £{exp (ja) Q t) + exp (-y'o) 0 f)} = cos w 0 t 
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Note that both lines have magnitude n which is in agreement with 
the general result a(— co) = a(w ) obtained above. 

It is worth while considering the spectrum of sine of. We wish to 
produce a real variation (i.e. along the x axis) which is sinusoidal. 
Again we produce the variation by adding two oppositely rotating 
vectors and they must be combined in such a way that the value of 
x at / = 0 is zero and, at first, x must increase as / increases. The vectors 
exp(/co/) and exp( -yco/) both have value 1 at / = 0. To make them 
cancel at this instant we multiply one by j and the other by - j , so 
that one lies along the direction of positive y and one along the 
direction of negative y. In order to make x increase with / at first the 
vector lying along the positive direction of y must rotate in the 
negative direction. Thus we require the sum - j exp ( yco/ ) 4- 7 ex p (- y'co/). 
If we specify unit amplitude then we must take half this sum. We note 
that this can be written 


{exp (jut) - exp (—jail)} 

which is the well-known formula for sin co/. Thus the spectrum in this 
case is pure imaginary and consists of two lines one of magnitude 
— 7 r at -f co and one of + it at — co. Note that this is the imaginary part 
of F(co) and that if we had plotted b( co) the signs would have been 
reversed, i.e. a positive value for positive frequency and a negative 
value for negative frequency since F(co) = a(aj)—jb(u)). Again we see 
that the spectrum being pure imaginary is in agreement with the 
general result deduced above for odd functions and also the magnitude 
of the two lines are equal and opposite, which is in agreement with the 
general result b(-w) = -£(co). 

We can combine the results for cos co/ and sin co/ to obtain the 
spectrum of the exponential function. Thus we can write symbolically 


fit) = sin co/, 


fit) — ysinco/, 


F(co) = a(aj) -jb(co) 


_ o Hines +7T at + a/1 

L — 7T at — w J 

+ 6(co) 

[ lines + 7r at + col 
— 7r at — co J 


F(co) = ja(a)) + b(aj) 
= y\0 + 
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fit) = cos co/, F(<jj) = a(cv) —jb(cv) 

[ lines + 7 r at + col , ^ 

+ „a,-»K 0 

The quantities in the square brackets in the last two relations are 
real and hence can be added (see section 3.8 on the superposition 
theorem). Thus 

fit) = cos co/ +7 sin co/, F(co) = [line+27rat+co]+y*.0 
= exp (yco/) 

Thus the complex spectrum of exp(y«>o0 is Pure real and consists 
of a single line of magnitude 2tt at + co 0 . This is readily verified for 

1 f® 1 

y-CO = _ j F{w) exp (Jut) dt = — [2tt exp (yo > 0 /)] = exp (yo » 0 /) 

Similarly the spectrum of exp (~j<u 0 t) is a single line of magnitude 
27r at co = — co 0 . 

3.7 The Shift Theorems 

We shall now deduce a relationship between the spectrum of a wave¬ 
form which has been shifted in time and the spectrum of the original 
waveform. Let F(co) be the Fourier transform of/(/) so that 

F(to) = /(/) exp (—j(ot)dt 

and let F_ r (co) be the Fourier transform of/(/- F), i.e. the original 
waveform delayed by T. Then 

F_ r (co) = J ^ fit — T) exp( —yco/) dt 

Putting t-T=t' 

F. jim) = /(f') exp { + T)}dt' 

= exp (-> 7 ’)J“ oo /(/') exp ( dt' 

Hence /' being a dummy variable we have 

F_ r (co) = exp (—yco7) F(co) 
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We can deduce formally a second shift theorem by a similar 
procedure, as follows. Let /(f) be the inverse transform of F(u) so 
that 

i r°° 

fit) = — I F(u) exp (Jcot) du 

and let f_ p (t) be the inverse transform of F(u -p) so that 

1 C 00 

f-pit) = — J F(u -p) exp (Jut) du 

Putting u —p = u 

f-pit) = J F(co') exp {j(t o' +p) t) du 

1 f 00 

= exp (Jpt) — F(u ') exp ( ju' t)du f 

2n J —oo 

Hence 

f-pit) = exp(y/>f)/(f) 

Expressing these two theorems in words: 

(i) If the spectrum of /(/) is F(u), then the spectrum of fit-T) 
is exp(-juT) F(u). 

(ii) If the signal represented by F(u) is /(f), then the signal repre¬ 
sented by F(u-p) is exp (Jpt) /(/). 

We may note at this stage that both of these theorems have been 
deduced in a formal mathematical fashion. The physical interpre¬ 
tation of the first shift theorem presents no difficulty since if/(/) is a 
real function of time then as already noted f(t—T) represents the 
same signal delayed in time by amount T. The second shift theorem is 
not quite so straightforward however. We saw in section 3.4 that for 
a real function/(/) the complex Fourier transform has the property 
F*(u) = F(— o>), i.e. the real part of F(u) is symmetrical (about the 
origin u = 0) and the imaginary part is anti-symmetrical. If now we 
substitute u—p for u then this relation will clearly no longer be true 
since in general the function F cannot have the same symmetry pro¬ 
perties about o) = p. Hence the signal represented by F(u - p) will not 
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be a real function of /. The truth of this is evident from the expression 
deduced for f- p (t) above. Although in general we shall not be 
concerned with complex functions of time we shall use the second 
shift theorem when we come to consider modulation theory. 

3.8 The Transform of a Derivative 

If we have a function f(t) whose Fourier transform is F(u) then we 
can readily obtain the transform of the derivative of/(f). We have 

i r°° 

fit) = — F(u)exp(jui)du 

27rJ_oo 

and if we let g(t) be the derivative of/(f) then 
d 1 d C 

g(t) = -/(/) = 2 ^J t \ F(a>)exp(jwt)d<o 

Reversing the order of differentiation and integration and carrying 
out the differentiation gives 

1 f 00 

g(t) = — {ju>F(u))} exp (ju)t) da) 

2.TT J — co 

We therefore see that g(t) and {j<*>Fi<*>)) form a Fourier pair. The 
Fourier transform of the derivative of/(f) is accordingly ju times the 
Fourier transform of /(f). The result is readily extended to higher 
order derivatives. 


3.9 The Superposition Theorem 

In this section we shall deduce an expression for the spectrum of the 
sum of two waveforms. Suppose the spectrum of fff) is F { ((o) and 
that of/ 2 (f) is F 2 iu). Then 

1 f 00 

flit) = — J F,(w) exp (Jut) du 


/ 2 (f) = — J F 2 (u) exp (jut) du 
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Therefore 

/.( 0+/ 2 (0 = ^J_“ 

+ 


F l (oj)cxp(J(x)t)du) 


i_ r 


F 2 (<o)exp(ja)t)d<D 


i r°° 

/iW +/ 2 (0 = — J (F^co) + F 2 (oj)} exp (Jcot) daj 


Hence the spectrum of {/,(/) +/ 2 (r)} is {Ffa) + F 2 (w)} or the 
spectrum of the sum of two waveforms is equal to the sum of their 
individual spectra. This is the superposition theorem. 


3.10 The Duality Theorems 

Let us consider a function of time/(/) and a function of frequency 
g(a)), and suppose/(r) and g(co) form a Fourier pair. Then we can 
write 

1 f 00 

/(*) = ~J g{oj) exp (Jwt) do) 

1 f 00 

• •f(-t) = — I g(a)) exp ( —yeti/) do) 

Hence by interchanging tu and t we obtain 
1 f 30 

/(- o>) = — I g(/) exp (-jcot) dt 

or 27r/(-c*>) = ^ #(f) exp ( -dt 

Comparing this with the general relation 

F((o) = f{t) exp ( -ja)t) dt 

we see that and g{t) are also a Fourier pair. 

Similarly we could proceed from the relation 

Efo>) = exp (-j<Dt)dt 

.*. g(-aj) = /(/) exp (yW) dt 
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and interchanging co and / 

foo 

g(~t) = I ^f((t)) exp (Jcot) dta 

so that — g(- /) and /(a>) also form a Fourier pair. 

2tt 

We can state the duality theorems in words as follows: 

(i) If the spectrum of f(t) is g(cj) then the spectrum of g(t) is 

2tt/(- cd). 

(ii) If the signal ofg(co)is/(/) then the signal of/(a>) is ^ -g(-t ). 

2tt 


3.11 The Multiplication Theorem 

We can obtain a formula for Fourier integrals corresponding to that 
derived in section 1.8 for Fourier series. Consider the integral of the 
product of two functions fi(t),f 2 (t). 

f f\(t)fi(t)dt = ^- f f\(i)dt f F 2 (w) exp C/co/) dm 

J — 00 ^‘ 7T J —oo J — oo 


-oo 

(*oo 


I /'oo /'oo 

= — I F 2 (to)du> /[(/) exp (Jwt) dt 

2 tt j —co J— oo 

i_ r 

" 277 J_ a 


Fj (—ai) F 2 (co) dm 


(w) F 2 (oj) dm 


Thus 


/'oo 

J-a 


mmdt = L 

Similarly it may be shown that 

J /l(0/2(0 d ‘ = ^ j ^ F l(“0 f 2*(") 
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Hence we have 


r mm* = ~ 

J —oo 27 T 


2 t t 



F* (cu) F 2 (co) dco 
F { (co) F*(co) dco 


3.12 Parseval’s Theorem 

As in the case of Fourier series we derive Parseval’s theorem by putting 
f\ =fi in the multiplication theorem. Hence we have 


or 


n = bS- p * (w) p{oj) dw 

/>»>* =— iFMpdo 


3.13 Energy Spectrum and Autocorrelation Function 

The total energy in a waveform is proportional to the integral of the 

square of its amplitude 

{/«»’<« 

which by Parseval’s theorem is equal to 

sjr„ 

The expression |F(a>)| 2 is referred to as the energy spectrum. Now 
by the multiplication theorem 

£A<M>* = £ F*(oi)F 2 iw)doi 

Let f 2 (t) =/,(/+r). Then by the first shift theorem the transform 
of/i(/+r) is given by exp (jojt) F,(cu) where F^oj) is the transform of 

/i«. 

J f\U)f\(t + r)dt = ^ J F*(w) exp (Jwr) Fj(oj) dco 


2tt 


/•oo 

| F, (w) | 2 exp (ycoT) dco 
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We define <f> x j(r) by the relation 

+r) dt 

The function <f> u (r) is known as the autocorrelation function. 
Hence we have 

i r°° 

*iW = 2^ I \F { (co)\ 2 exQ<J(OT)d(o 

Thus 4>\\ir) and |F,(cl»)| 2 form a Fourier pair, and the energy 
spectrum is given by the Fourier transform of the autocorrelation 
function. 

The autocorrelation function is of great importance in the analysis 
of random waveforms and also in probability theory. However a 
discussion of applications in these subjects is beyond the scope of the 
present work. 


3.14 The Convolution Integral 

Consider the Fourier transform of the product of two functions fi(t) 
and f 2 (t). If the Fourier transforms of f Y (t) and f 2 (t) are F x {(o) and 
F 2 (w) then the transform of the product, F(co) may be obtained as 
follows. 

F(co) = j mm exp ( -]<ot) dt 





/ 2 (0exp(-yo>/)<// Fi(p) exp (jpt) dp 


r oo 

Fi(p) dp f 2 (t ) exp { —j(a> —p)t)dt 

J — 00 

Flip) F 2 ((i>—p) dp 


and so we have the result 


/•« 

J- < 


flit)fiit) exp ( -jwt) dt 


j_ r°° 

_ 2 w J — 0 


Fi(p)F 2 ioi-p)dp 
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The integral on the right hand side is known as the convolution 
integral. Thus the spectrum of the product of two functions is the 
convolution of the individual spectra of the two functions. The 
convolution operation is sometimes written symbolically 

F{<*>) = ~~ F 1 (<o)*F 2 ( oj) 

17T 

We can also derive the dual of this theorem by taking the inverse 
transform of the product of two spectra. 



COO 

Fl( co) F 2 (to) exp C/cot) dco 

— co 

-sj 

1* CO 

F 2 (co) exp (ycor) dej 
' — 00 « 

C 00 

1 f Y {t)exp{-]o}t)dt 

II 

N> 1 _ 

3 | — 

«_ 

p oo Coo 

fl(t)dt F 2 {oj) exp {y'co(T - 0} dm 

— co J — CO 

= r mh(j-t)dt 

1 m 


and so 



J Fi(co) F 2 (w) exp (yW) da> = J 

’ mf 2 {r-t)dt 
— 00 


= /i(0*/ 2 « 


The signal of the product of two spectra is thus the convolution of 
the individual signals of the two spectra. 
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CHAPTER IV 


Analysis of Transients 


4.1 Introduction 

Just as Fourier series can be used in the analysis of periodic waveforms 
as seen in Chapter II so Fourier integrals can be used in the analysis 
of transient waveforms and in this chapter we shall consider some 
examples. 

4.2 The Rectangular Pulse 

Consider the rectangular pulse shown in Fig. 4.1(a) defined by the 
relations 

/(')=■ 2 r ~T<t<+T 
= 0 elsewhere 


Note that this definition results in a pulse of unit area. We obtain 
the spectrum of the pulse by taking the Fourier transform 


F(oj) 


Coo 

J- 0 


f(t)exp(-jwt)dt 


.±r 

IT}., 


exp(-yco/)<// 


J_(_zl) 

2T jej 


[exp (-yco/)] r r 


or 


F(<o) = 


{exp C/wT) - exp ( -jwT)} 

sine oT 
cjT 


This expression will be zero whenever sino>T= 0, i.e. for to = nn/T 
except at the origin when the denominator becomes zero also. When 
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a>T is small sincoT is approximately equal to ojT and the value of 
F(cj) at (o = 0 becomes unity. As w becomes larger sin ojT continues 
to oscillate between +1 and —1 whereas cjT the denominator 



Fig. 4.1 

(a) The rectangular pulse and ( b ) its spectrum 


increases continuously. Hence as cj increases the magnitude of the 
oscillations decrease and finally become zero. A sketch of the spectrum 
is shown in Fig. 4.1 (b). 


4.3 The Rectangular Frequency Distribution 

We shall consider now the analogous case of a rectangular frequency 
distribution (Fig. 4.2) defined by 



Fig. 4.2 

(a) The rectangular frequency distribution and ( b ) its signal 


We have 
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1 f°° 

AO = 2n) ex P dw 
1 1 C°“ 

= — — exp (jut) dm 

2w2w c J_ m< 

= ^4— {exp CK 0 - exp ( -j<o c t)j 


1 sin oj c t 
2 tt a) c t 


We could have obtained this result from the duality theorem and 
the results of the last section. The duality theorem states that if the 
signal of g(cu) is /(f) then the signal of f(cS) is g(-f)/27r. Thus we 
define g(a>) by the relation 


g(<o) = 


sin cuT 
ojT 


and 



-r<f < +t 

elsewhere 


Thus/(f) is the signal of g(co). The function f(<o) would then be 
defined by 


/(<*) = 


2co c 


— 0) c <UJ<W c 


= 0 elsewhere 


where we have replaced the particular value of time, T, by the particu¬ 
lar value of frequency o> c . Hence by the duality theorem the signal of 
this spectrum is given by 

1 1 sin(-o> c f) _ 1 sin <o c t 

2tt^ 2tt (— o > c f) 2tt w c t 

where again Tis replaced by a> c . 
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4.4 Reciprocal Spreading 

We have seen in the discussion above that the first zero in the fre¬ 
quency spectrum of the signal confined to the region — T< t < T in 
time occurs at a frequency = 7 t/T. Since the majority of the signal 
energy will be contained in frequencies less than this we can consider 
the bandwidth of the signal to be of this order, i.e. we may write 
/V 1 /2 T. The time A t occupied by the signal is 2T and hence we have 

Af.At ~ 1 

Hence as the region in time occupied by the signal becomes smaller 
so the region it occupies in frequency becomes larger. A similar 
result is obtained by considering the signal confined to the region of 
frequencies — a) c < w < a> c , i.e. frequencies less than w c so that 
Af= wJItt. Here the signal energy lies mainly in the time interval 
-7r/u) c <t<7T/cj c so that At = 2 ttIo) c again giving Af.At~ 1 . A 
signal can be represented either as a function of time or as a function 
of frequency (i.e. its spectrum) and as it is compressed in one repre¬ 
sentation so it expands in the other. This phenomenon is sometimes 
referred to as reciprocal spreading. 


4.5 The Unit Impulse 

This is an important function and can be developed from the rect¬ 
angular pulse of section 4.2. This was defined in such a way that its 
area was unity. The unit impulse is obtained by taking the limit 
T-+ 0. The amplitude of the function then becomes infinite at the 
origin but the area remains unity. The spectrum retains the value one 
at the origin w = 0 and the first zero on each side of the origin move 
out to infinity. The spectrum then becomes 


F(co) = lim 
t —> 0 


sincuT 

ojT 


= 1 


i.e. the amplitude of the spectrum is independent of w and becomes a 
straight line parallel to the o> axis at amplitude one. The energy of 
the unit impulse is infinite. 

By the duality theorem or from section 4.3 we deduce that a 
frequency impulse at the origin gives a signal of constant height 1 '2tt, 
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i.e. a frequency impulse at the origin represents a d.c. signal, which is 
obviously correct. 

Now applying the second shift theorem a frequency impulse 
occurring at frequency a>=p will represent a signal exp (//?/)/(/). 
This, as was explained in the discussion on the shift theorem is not a 
real function of time. We can overcome this difficulty however as we 
saw in section 3.6 by adding a second impulse at frequency w = — p 
which by the shift theorem gives a signal exp (-jpt)f(t). Since/(/) was 
a constant ( 1 /2tt) the resultant signal is (1 /7r) cos pt, by the superposition 
theorem. Hence a sinusoidal type of variation gives rise to frequency 
impulses or spectral lines, a result we are already familiar with. Note 
that for the signal cos pt the magnitude of the spectral lines will be 7 r, 
for 

1 

fit) = — I F(w) exp O’o>0 dw 

= ex P (Jpt) +7 T exp ( -jpt)} = cos pt 

27 T 

as required (see also section 3.6). 

4.6 Finite Train of Pulses 

From section 4.2 we have that the spectrum of a pulse of unit area 
and duration 2t at the origin is sin(an-)/(coT). Thus a pulse occurring 


fit) 

rc 


0 t _ > T-T T T+T 

Fio. 4.3 

A rectangular pulse occurring 
at time T 


at time t — T as shown in Fig. 4.3 will have a spectrum given by 
exp0‘coT). sin (o>t)/(ojt), by the shift theorem. Now consider a train 
of N pulses occurring at 

/ = -kT 9 -(k-\)T,...-T 9 0,T,...(k-\)T,kT 
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where N= 2Jc + 1. Then by the superposition theorem the spectrum 
is given by 

, v sin cot. ^ _ 

F( co) =-[1 4 - exp (jwT) 4- exp (-ja>T) 4 -... 

COT 

4 - exp (JkwT) 4- exp ( -jkwT)] 

=-[— 1 4 ~ 2 {l 4 * cos coT 4 - cos 2 coT +... 4 * cos kwT^\ 


The series in square brackets may be summed by taking the real 
part of the corresponding exponential series, i.e. 

-14- 2 Re [14- exp (JwT) 4-... 4- exp (jkwT)\ 

f 1 — exp {j(k + 1 ) 

L 1 -expO'^r) J 

|~ {1 ~ exp M + 1 ) a>r}{ 1 - exp (-ycoDj l 

2(1 — cos wT) J 

1 - cos (wT) 4- cos (kcoT) - cos (k 4 -1) wT 


= — 14- 2 Re 


= — 1 4 * 2 Re 


= -1 + - 


1- cos wT 


cos (, kcoT ) — cos (k 4 -1) c oT 
1 — cos wT 


2 sin i(2k 4 -1) co7\ sin \<jjT 
2 sin 2 \u)T 


sin \NwT 
sin \wT 


Hence we have 


F(co) = 


sin cor sin J NwT 
cot sin }coT 


which is the spectrum of the finite train of pulses. 

Let us now consider the behaviour as N-+ oo. We can write 


/ro¬ 



sin cot sin \NwT 
cot sin i<oT 


exp ( jwt ) c/co 
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Using again the theorem J« <£(£)sina£c/£-*0 as a->oo provided 
<£(£) is finite in the region a < £ < b we see that the only contribution 
to the integral will be where the function 


<£(co) = 


sin cot 


1 


cot sin \wT 


exp (Jwt) 


is not finite, i.e. where 


sin \wT = 0 


This occurs when 


\wT = wr or co 


2tm 

~T 


The contribution from one such point will be 



f(2 ™IT)+* sin (IttutIT) [ limAsm J NwT J 

-exp I /- 1 -c/co 

'(2m/d-« (2irnr/T) \ T J sin \a>T 


where ui has been replaced by the value 2-irnlT in those functions 
which are finite and slowly varying in the region of integration. 




J_ sin (2 nnrlT) / 2nnt \ f <2 ™ /r)+a sinjAfar 
2ir (2nnrlT) y T / J^n/D-a sin i<oT 


Now let INcoT = £. For large N we have 


/•(2 oti/d+« sinJAilwT’ j*® sin | 2 f“ sin| 2 

Ja™/D-a sir>ia>r ~ J-* sia((/N) NT ^ ~ J - n W)NT' 



sinf 2 rr 

r 


:-m = 


1 sin (27 rnr/T) 
- c 

27r (JjttvtIT) 

1 sin (limr/T) 
2tt tit 


/ 27rnt\2ir 

[ J ~i -)t 

( 27rnt\ 

M 


c 
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Hence/(r) is given by 

1 sin (277/ir/r) 


/«>- 2 


2 rr 


/IT 


1 [27T 'SpSL 

“ 2 ^|T + 2 ” 

,/«. i + i I 


sin (Ittht/T) \ 


m 


l lirnA 


sin (Inrir/T) 27/ 

-cos —nt 

tit T 


This result could of course have been obtained more easily by the 
Fourier series method 




1 2tt 

— cos — ntdt 


Hence as before 


1 sin ( lirnr/T) 
7r nr 




sin C Ittht/T) 
m 


2tt 

cos — nt 
T 


4.7 The Gaussian Function 
Consider the simple form of the Gaussian function 
/(/) = exp (-Jr 2 ) 

The Fourier transform is given by 

F(co) = J ^ exp(- it 2 )exp(-jajt)dt 
= exp {—i(f 2 + 2j(ot)} dt 

= exp ( -\a> 2 ) exp {-£(/* + Vut~ 

- exp(- W) exp { - i(t +» 2 } dt 

= V2exp(- W) J* w exp { - i(t +joj) 2 } d{(t +»/V2} 
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But 

J 00 ^ exp (-x*)dx = y/jr 
F(a>) = V(2*0exp(-i<o 2 ) 

Thus the Gaussian function has the same form, viz. exp(—ix 2 ) 
in both the time and frequency representations (cf. the infinite train 
of narrow pulses considered in section 2.6) and is said to be self- 
reciprocal. 

It is interesting to note that the Gaussian function is one of a set 
known as the Hermite functions. These are given by the relation 

H n (f) = ^ exp (-J/ 2 ) 

Thus the Gaussian function is Hq(J). All the Hermite functions have 
the property of self-reciprocity. 

4.8 The Unit Step Function 

The unit step function is defined by the relations 

/»- 1 />o 
= 0 / < 0 

and is shown in Fig. 4.4. This is another important function which 
finds many applications in engineering, e.g. in circuit theory. We note 

rtt) 

i —- ■■■■— 


o t-> 

Fig. 4.4 

The unit step function 

that as defined the unit step function does not satisfy Dirichlet’s 
conditions since|/(/)| dt is not finite. To overcome this difficulty 
we define a new function 

/(/) = 1 0 < t<T 
= 0 elsewhere 
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Dirichlet’s conditions are now satisfied and we can proceed with 
this function. The spectrum is given by 

F(u>) = /(/) exp (-jwt) dt = J Q r exp (->/)<* 


= t- {l — exp ( —jtoT)} 


Hence we have 


i r°° i 

fit) = — J j—{l - exp i-jcoT)} exp (/a>/) da> 

We can make this function approximate the unit step function by 
making T large and under these conditions 


J 00 \ r*> 

. J__ 


cxp(-jw T) 
jw 


dw 


which is a standard integral and has value -v for T> 0. 

exp (jwt) 


m 


-iC 


j(X) 


-daj 


This result is readily checked since the integral is similarly a 
standard integral and 


J. 


exp (jwt) 


-oo J<*> 


d(i) = — 7T / < 0 


Hence 


fit) = 1 
= 0 


= +7T t > 0 
t > 0 


as required. 


/ < 0 
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From the formula 

/ (0 = i + ^ J exp (Jwt) du> 

we see that there is a d.c. term of magnitude i which has a spectrum 
consisting of an impulse of magnitude n situated at = 0 (see section 
4.5). The spectrum of the second term in the expression is clearly 
(1 Ijcj) and so from the superposition theorem the complete spectrum 
of the unit step function is an impulse of magnitude tt at co=0 plus 
the function F(aj) = (l//a>). Hence the magnitude of the component 
frequencies in the unit step function is inversely proportional to the 
frequency. 

4.9 The Exponential Decay 

This is another function which frequently arises in engineering 
practice. It is defined by 

/(/) = exp(-a/) / > 0 

= 0 / < 0 

A sketch of the function is shown in Fig. 4.5. The determination of 
the spectrum is quite straightforward. 



An exponential decay 
F(w) = /(0 exp dt 

= J” exp (—at) exp ( —jwt) dt = J" exp {-(«+» i]dt 


-1 
a +ja) 


[exp {-(a+»*}]£ = 


1 


a+yuj 
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Thus 


F(<o) = 


1 

a +jco 


|F(«>)| = /=■*(«>)] = 


J[ 


1 


(a +ycu) (a 



A |F(a>)| = 


1 

V(a 2 + a> 2 ) 


4.10 The Exponential Rise 


This function is related to the previous one and they often occur 
together in different parts of the same system (e.g. the transient 
response of an RC circuit). It is sketched in Fig. 4.6 and defined by 


/(/) = l-exp(-af) / > 0 

= 0 t < 0 



An exponentially rising 
function 


This function, like the unit, step does not satisfy the condition 
|-oo |/(01 dt finite and its spectrum is most readily obtained by means 
of the superposition theorem. We can consider this function as a unit 
step minus the function exp(-ocf). The spectrum of the unit step is 
given by F(co) = 1/y‘cu and the formula for/(/) is 

1 f 00 1 

fit) = i + — I r-zxvUwOdu) 

Hence for the exponential rise the spectrum is given by 


joj a + jco ycu(a + jw) 
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and 


and 


l F(<0)| = a>\/ (a 2 + co 2 ) 

m J_ p _a_ 

27rJ- m j<o(<x+j< 


f(0 


exp (jcot) dto 


4.11 The Cosine Wave of Finite Duration 

This is obviously another case which arises in practice since it repre¬ 
sents a co-sinusoidal oscillation which is switched on at some point 
in time and then off again at a later time. The r.f. pulse commonly 
used in radar engineering is just such a waveform. A sketch of the 
function is shown in Fig. 4.7. This function is defined by 

/(/) = cos pt —\T<t<\T 

= 0 elsewhere 



A cosine wave of finite duration 


We obtain the spectrum as 

r® rr /2 

F(ao>) = f(t)exp(—jcot)dt = cos pt exp (-jwt)dt 

J—oo J — 7/2 

Whilst this is a standard integral and the integration can quite 
readily be carried through it is a little lengthy to do. A simpler integral 
results if we make use of the properties of F(co) which we have 
deduced. The function f(t) in this case is even and therefore its Fourier 
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transform is a pure real function, i.e. it is equal to a(co) the cosine 
transform (section 3.4). Hence 

F(co) = a( co) = cos cot dt 

J r /2 

_ TI2 cos co/cos/?/cfr 

rT/2 

= i J _ r/2 {cos (co -p) i +cos (co +p ) /}dt 


= i[— — sin (co-/?)/+— sin (co +/?)/] 
[co-/? co+/> J 

1 • , 1 . , T 

--sin (co-/?)-H -sin (co +/?) - 

co-/? 2 w+p 2 


Til 

-TJ2 


which may be written 


F( , _ Th\nl(w-p)T sinKw+p)rl 

2[ «o >-p)T + Kw+p)T J 

The number of cycles of the waveform in the pulse, i.e. the number 
of cycles of angular frequency p occurring in time Fis given by pT/lir. 
If there are a large number of cycles in the pulse, i.e. pT is large, then 
the second term in the expression for F(co) will be small for the 
positive range of co. The curve for F(co) then becomes the curve 
sin*/* centred at co=/?. When however there are only a small 
number of cycles in the waveform either because p is small, i.e. a low 
frequency, or because T is small, i.e. a short pulse, then the second 
term makes an important contribution to F(co) and the curve changes 
its shape. 

Consider for example the case shown in Fig. 4.8 (a) where there is 
only one half cycle in the waveform. In this case pT= ?r. Hence we 
have 


17/ \ 1 /tTCO 77" \ 

F(co) = -sm --- - 

co-/? \2 p 2) 



• EY \ ~ 2 P n* 0 

.. r (co) = -z- =cos-- 

cx) z —p 2 2 p 
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and on writing co//? = x this becomes 


F (< o ) = 2 - 3 -. COS tirx 

Pil-X 2 ) 



(a) A cosine wave of half period duration and (6) its spectrum 

We can sketch this curve as follows. For co = 0, i.e. x = 0, we have 
F(0) = 2Ip. At the point co =/?, or x = 1, both numerator and de¬ 
nominator become zero and we evaluate the expression by differ¬ 
entiating numerator and denominator and then putting x = 1. Thus 



The expression for F(co) has zeros where cos(£7tjc) = 0, i.e. where 
£77 x = £(2/i+l)7r, or co = (2 n+ 1 )p (except the point co = p which we 
have already dealt with. When x becomes large cos(£7tx) oscillates 
between +1 and -1, and (1 -x 2 ) increases in magnitude continu¬ 
ously. Hence the magnitude of the oscillations of the expression for 
F(co) decreases continuously as co increases. The curve is sketched in 
Fig. 4.8(6). 

It is instructive to examine this change of shape a little more closely 
We derived the expression for F(co) 

= Fr sin£(co-/?)F sin£(co+/?)F l 

{W) 2[ 4 (co-/?)F 4 (*+p)T J 

Thus the expression consists of the sum of two curves of the form 
sinx/x, one centred at co = +/? and one at co = -p . We might have 
expected this from our investigation of the impulse in section 4.5. 
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We saw there that an impulse at the origin represents a constant 
d.c. from — co to + oo. If however instead of a constant value, i.e. a 
frequency of zero for the time function we have a frequency p then the 
spectrum splits into two impulses of half the magnitude spaced 
symmetrically about the origin of w (i.e. at oj = ±p). We know from 
section 4.2 that if the finite value of the d.c. is limited to a finite 
interval of time instead of extending to infinity then the spectrum 
takes the form of the curve sin jc/jc centred on the origin. Hence if we 
have a frequency p in this region instead of zero (d.c.) we might have 
expected the spectrum to split into two curves of the same shape 
(i.e. sinx/x) centred at cj = ±p. This can readily be verified by the 
shift theorem. Let f(t) be defined by 


f(t) = 1 -\T<t<\T 

= 0 elsewhere 


Then from section 4.2 


F(a>) = T 


sin \o)T 
l<oT 


If now we shift this curve to w = p and take half the magnitude then 
we have 


F(w—p) 


T s\n\(u)-p)T 
2 ~l(c»-p)T 


By the second shift theorem the signal of this spectrum is given by 
i exp (//?/)./(/). Shifting the curve to <u = — p gives 


Fico+p) = r[!l£i^MZl 

2[ b((o+p)T J 


which by the shift theorem has signal Jexp (-jpt )./(/). Hence by the 
superposition theorem the signal due to a spectrum 


F(<o) = sin j(cu-f p) T ~j 

21 i(o>~p)T J («>+p)T J 

is given by 

i{exp (Jpt) + exr (-jpt)}f(t) 

= cospt.f(t) = cos pt - iT < t < IT 


= 0 elsewhere 
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We can use these ideas to explain the change of shape of the curve 
noted above. When there are a large number of cycles in the waveform 
pTis large and the two branches of the curve of F(cj) are independent. 
The one centred at o> = -/> becomes negligibly small over the region 
where the one centred at <d = +p is important as shown in Fig. 4.9(a). 




Fig. 4.9 

Change of spectrum as the duration of the cosine wave varies: (a) the positive 
and negative parts are isolated from each other; ( b) some overlap occurs; 
and (c) the case for half period duration 

If the time Tis reduced keeping p constant then the positions of the 
first minima move away from the point at which the curve is centred 
and the curves are stretched out in the / direction. If p is reduced 
keeping T constant then the two branches of the curve move towards 
each other keeping the same shape. In both cases a point is eventually 
reached when the two curves begin to overlap and the curve centred 
at a) = —p still has an appreciable value over part of the region where 
the curve centred at cj = +p is important. This condition is shown in 
Fig. 4.9(6). Since the curve for F(w) is the sum of the two branches 
then at the point where the overlap becomes important the curve 
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begins to change shape. The condition pT= tt, the case considered 
above, is shown in Fig. 4.9(c). In this case the two curves combine 
in such a way as to give a resultant curve with one main peak at the 
origin. From the figure we see that the value of this peak is given by 

[ 2 Tsin*1 _ IT _ 2 

^ X J.v=$ir 7T P 
since for this case pT— it. 

Returning now to the case shown in Fig. 4.9(a) it can be seen that 
the majority of the energy of the signal is contained in a range of 
frequencies Afoi the order 1 /T. The time duration of the signal we 
considered was T. Hence again we have the relation obtained in 
section 4.4, viz. Af.At~ 1 . If we wish to define the frequency of a 
signal accurately, i.e. to make Af small, then At must be large. 
Conversely, to make At small we require a large bandwidth. Since a 
signal of finite duration must have a finite bandwidth its frequency 
will be to some extent uncertain, and the relation Af.At~ 1 is 
sometimes referred to as the uncertainty relation. 


CHAPTER V 


Application to Circuit Analysis 

5.1 Introduction 

In this chapter we shall apply the ideas developed earlier in the book 
to some of the problems which arise in dealing with electrical circuits. 

5.2 Filters 

We shall consider in this section only linear filters, i.e. those filters 
whose output is proportional to the input. The factor of proportion¬ 
ality however is not the same at all frequencies, and this gives rise to 
the filtering action. The filter may be specified in terms of its frequency 
characteristic or transfer function as it is sometimes termed. Thus if we 
represent the transfer function by Z T (cu) then the output at frequency 
<x) 0 is Z t (cj 0 ) times the input. Hence if the input is represented as a 
function of frequency, i.e. as a spectrum F^to), then the output 
represented as a function of frequency is 

F» = Zj{w)Ffa) 

If, then, we have a waveform /,(/) as the input to a filter we can 
obtain the output by the use of Fourier transforms. We can write each 
of the complex functions in the above equation in terms of amplitude 
and phase. Hence we have 

R 0 (w)cxp {j0 o (a))} = A(w) exp {j<f>(o>)} F*(o>) exp {jOfa)} 

where A{oS) is the amplitude characteristic of the filter and </>(co) its 
phase characteristic. 

It is convenient to consider the general case where a) is allowed to 
take negative as well as positive values and we shall now deduce some 
of the properties of Z T (w). We know that 

F(-a>) = F*(co) 

and so it follows immediately that 

Zj>(—ou) = Z*(w) 
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and hence A(-oj) = A(w) 

and </>i-<») = 

The amplitude characteristic is thus an even or symmetric function 
of w and the phase characteristic is odd or anti-symmetric. 

5.3 Distortionless Transmission 

If a waveform is to be transmitted without distortion then the only 
operations which can be carried out upon it are to multiply it by a 
constant and to shift it in time. These operations correspond to 
amplification or attenuation, and delay. If we multiply a time function 
fit) by a constant A then all the frequency components are also 
multiplied by A. If the function is shifted in time, i.e. it becomes 
f(t - T) f then the frequency representation is multiplied by expi-jwT), 
by the shift theorem. Thus the general characteristic of a filter which 
will transmit a waveform without distortion is 


Z T i<j>) = Aexpi-jwT) 

This is readily verified, for if the input wave is given by fit) then 


and 
so that 


1 f°° 

fit) = — I Ffiw) exp ijwt) dw 

F 0 iw) = Z T {a>) Fi(a)) = Aexp(-jwT)F { iw) 


fait) 


F 0 (o>) exp ijo)t) do) 


— f 00 
2tT J_a 

i r°° 

" 2 ^ I {/I exp ( -jo)T) Ffw)} exp ijwt) dw 


or fait) = Afit-T) 

We can illustrate this by considering a transmission line. If Vix,t) 
is the voltage at a point distant x from the input at time / then it may 
be shown that for a sinusoidal input of frequency w 0 


Vix,t) = tT expj(w 0 t—Px) 
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where V is the amplitude of the voltage at the input and P is the 
propagation constant. The input voltage is given by K,(/) = F(0,/) 
and the output voltage V Q it) = ViX y t), where X is the length of the 
transmission line. Thus we have 

v,(f) = ^expOoO 

and 

Voit) = 'T'expjiw 0 t-PX) = Viit) exp i-JPX) 

Thus at the frequency w 0 the output voltage is exp i-jPX) times 
the input where P the propagation constant will be a function of 
frequency. The magnitude of the frequency components at the output 
is obtained by multiplying those at the input by exp{-jPiw) X}. This 
function then is the transfer impedance of the transmission line. It may 
be shown that the function Piw) is given by 

P{w) = ±\/{— iG+jwC)iR+jo)L)} 

where G is the shunt conductance, C is the capacitance, R is the series 
resistance and L is the inductance per unit length of the transmission 
line. Hence we have 

Z T iw) = exp [ ±jXy /{—(G +jwC) iR +jwL)}] 

For distortionless transmission the transfer function must be of 
the form A exp (- jwT) and we can show that this form is obtained by 
making 

(G+jwC) = [M 2 iR+jwL ) 

where p is a constant. 

If this relation is satisfied then 

Zjiw) = exp [±jX\/{-n 2 iR+jwL) 2 }\ 

= exp [ ± pXiR +jwL)\ 

Z T iw) = exp (± fxXR) exp ( ±jfiXwL) 

We note from above that /z has dimensions [R]~ l and that R and 
L are measured per unit length. Thus fiXR is dimensionless and 
constant. Considering for example the negative signs we can write 
exp i-fiXR) = A. Similarly it can be seen that \lXwL is dimensionless 
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and can therefore be written a >T where Tis the time given by fiXL and 
is constant. Hence 

Z t (cj) = Aexp(-jwT) 

which is of the form required. Thus a transmission line will transmit 
signals without distortion if its parameters are related by the formula 

(G+jaiC) = fi 2 (R + jwL) 
a well-known result in transmission line theory. 

5.4 Passage of a Step Function through an Ideal Filter 
If a filter does not have a transfer characteristic of the form developed 
in the last section then, in general, distortion will result. For example, 
if A is a function of frequency so that some frequencies are attenuated 
more than others then the function of time appearing at the output 
will be different from that at the input. We shall investigate this 

AM 


-0) c 0 D c 

Fio. 5.1 

A perfect filter characteristic 

phenomenon by considering the passage of a step function through 
an ideal filter. We suppose that the filter passes all frequencies up to 
frequency a> c without modification and completely stops all those 
above co c . The transfer characteristic is therefore as shown in Fig. 5.1. 
This may be written 

A(w) = 1 — co c < co < co c 

= 0 elsewhere 

Consider the symmetrical step function given by 
fi(0 = t < 0 

= + J / > 0 
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Then from section 4.8 


Hence 


FiM = - 

JOJ 


F 0 (w) = Z r (co) Fi(aj) = — - co c < co < co c 

jw 

= 0 elsewhere 


so that 


m = 


JO) 

cos co/ 


_l r° i 

27r J_o, c yc< 

_l r 

2tt 

i r *- 

* Jo 


J<*> 
sin co/ 


exp (/co/) do) 
doj + 


1 r Wfl /sinco/ 

J-o, # j<*> 


do) 


do) 


since the first integral has an odd integrand in co and the second has 
an even one. Putting k = o)t 

r sin k 


m 


i ri 
* Jo 


-dk 


= -Sifat) 


where Si(o> c t) is the sine integral which we discussed in section 2.8. 
If we had considered the unit step function then we should have had 
the relation (see section 4.8) 


m = 



exp (Jot) 
/co 


do) 


Hence the output function f Q (t) would be as shown in Fig. 5.2. 

It is seen that the output voltage overshoots the final steady value 
and oscillates about it, the amplitude of the oscillations decreasing 
with time. As co c increases the value of the overshoot is not altered but 
the time at which the maximum occurs is smaller, the whole curve 
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being compressed towards the origin but remaining unaltered in the 
y direction (cf. Gibbs’ phenomenon, section 2.8). 

It is also seen that there is in addition to the overshoot what has 
been termed an anticipatory transient. Although the step function 
does not occur until the time / = 0 there are oscillations in the output 
waveform before this instant. We therefore have an output with no in¬ 
put in the region t < 0—the output anticipates the input. In the above 
analysis we neglected the phase characteristic of the filter. This as we 



The output of a perfect filter 
for a step function input at 
t =0 


have seen introduces delay. Hence taking this into account the whole 
curve of Fig. 5.2 should be shifted to the right. However the curve of 
Fig. 5.2 will not become zero (except at discrete points) until / becomes 
infinite and hence there is always a small anticipatory transient for a 
finite delay. A filter with the characteristics assumed is not physically 
realizable. 

5.5 Passage of a Step Function through a CR Circuit 
We have developed a method of calculating the output of a filter by 
first obtaining the spectrum of the output and then taking the inverse 
Fourier transform. In the last section we obtained an output which 
extended to infinity in both positive and negative directions. Its 
spectrum however was confined to a finite region of the frequency 
spectrum by virtue of the action of the filter assumed. This is an 
example of the phenomenon of reciprocal spreading which we dis¬ 
cussed in section 4.4. If the function is limited in one representation 
it spreads in the other. Thus to obtain a function of time which is 
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confined to a limited interval of time, say t > 0 , corresponding to a 
physical case, then the frequency representation must extend to 
infinity and can have no sharp cut-off of the type assumed in the last 
section. We have already calculated the Fourier transforms of several 
functions which have a discontinuity at the origin, i.e. functions 
which start at some definite point in time, a typical one being the 
exponential decay of section 4.9 which has spectrum l/(a-h/w). It 
is interesting to work through the inverse process to illustrate how 
spectra can be formed by filters which lead to time functions with a 
discontinuity at the origin. 

Consider for example the CR circuit shown in Fig. 5.3. The transfer 
function of this circuit is given by 

^ R jo) 

= R + (\ljmC) = a+jw 

where a = \/CR. 

0- 1 | -- 0 

Vi(t) C V,(t) 

Q 1 ■ - o 

Fig. 5.3 
The CR circuit 


If the input V t (t) is a step function then the output spectrum is given 
by 


1 ja) 1 

F 0 (w) = -— J — =-- 

ya> a +yca a 4 ja> 

We already know from the results of section 4.9 that the signal of 
this spectrum is given by 

fit) = exp (-at) t > 0 
= 0 t < 0 


and the output is therefore an exponential decay, a well-known result 
in circuit theory. It is interesting however to see how this discontinuity 
arises from the inverse transform. The output waveform is given by 




exp (Joit) 

- day 

oc+jo) 
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This integral can be evaluated by the contour integration method. 
Consider the integral 


f(z) = 


c ^ 

Jc *+JZ 


where C is a contour in the complex z plane. This type of integral is 
evaluated in terms of the residues at the poles, and in this case there 
is one pole at z = jot, which has residue exp (- a/). The contour Cmust 
be chosen so that the integral converges. In this case the infinite 
semi-circles on the real axis are appropriate, and which of these two 
is used depends upon the parameter /. If t is negative then for the 
integral to converge on the semi-circle at infinity the imaginary part 
of z must be negative and the semi-circle in the lower half plane must 
be used. This does not enclose any singularities of the function and 
hence the integral is zero. For t positive, then for the integral to con¬ 
verge on the semi-circle at infinity the imaginary part of z must also 
be positive and hence the semi-circle in the upper half plane is to be 
used. This encloses the pole at z— jot and the value of the integral is 
2tt times the residue at the pole, i.e. 27rexp(—a/). Thus the integral 
gives rise to a function of t which is zero for negative t and has a 
discontinuity at the origin. On completing the calculation the result is 
obtained 

f Q (t) = exp(-a/) />0 
= 0 t< 0 


in agreement with that obtained earlier. (For the details of carrying 
out this type of integration see, for example, Phillips (1943).) 

The transient response of circuits (the output when the input is a 
step function or an impulse function) is not usually calculated by 
means of Fourier transforms but by the use of another integral 
transform, the Laplace transform. In this case the kernel 


AT(p,/) = exp (-/?/) 

and the limits are from 0 to co. Thus the Laplace transform is defined 
by 

-^(/O = f(t) exp (-pt)dt 
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where p is in general complex. The exponential function now has a 
real part as well as an imaginary part and problems of convergence 
are less likely to be met with. This has proved a powerful tool in 
many applications and there are many books on the subject (see, for 
example, Jaeger (1961)). 


5.6 Initial Rate of Rise and Bandwidth 

We noted in section 5.4 that as the value of the cut-off frequency co c 
of the filter increases the curve for the time function is compressed 
towards the origin. Thus the rate of rise at the change-over increases. 
We can relate the initial rate of rise to the bandwidth co c as follows. 
We have for the output function 


m = 


1 f"«'sinA: 

—r dk 

* Jo k 


Now it is shown in the appendix that 


d_ 

dt 



g(x)dx = ag{at) 


Hence using this formula 


d a> c sin o) c t 

J, /oW 

dt 7T OJ c t 


This function is a maximum at / = 0 and hence the initial rate of 
rise is the maximum rate of rise also, and is given by 




7T 


We see that the rate of rise is proportional to the bandwidth and 
so the more rapid the rate of change in the signals which are to be 
transmitted, the wider the bandwidth of the filters must be. 

We can also calculate quite simply the value of the initial rate of rise 
for the CR circuit. Here we have 


foif) = l-exp(-a/) 
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The filter characteristic is given by 


Zj-(co) — 


1 

1 + jcoCR 


In this case the definition of bandwidth is somewhat arbitrary. If 
we take the ‘ 6 db point’, i.e. where \Z t (oj)\ = J then <d c CR = \/3 or 
<o c = VV(CR) 



CR V3 


This is a similar result to that obtained above for the ideal filter. 
Numerical agreement could not have been expected owing to the 
arbitrary definition of bandwidth co c in the second case. 

The foregoing discussion shows that the bandwidth required in a 
transmission system is closely related to the rise times it must transmit. 
The above analysis is readily extended to the case of pulses and here 
the rise times will determine how accurately the pulse shapes are 
preserved and also set the limit on the ability of the system to resolve 
consecutive pulses. It may be shown that two pulses separated by a 
time At will just be resolved by a filter with a bandwidth J/if J/and 
A t are related by the uncertainty relation Af.At~ 1 (see also sections 
4.4 and 4.1 1 ). Alternatively, the bandwidth required in the system can 
be obtained from a consideration of the spectrum of the pulses. A 
detailed discussion of these considerations will be found in Goldman 
(1948). 


5.7 Response of a Circuit to an Input Function 

Consider first the case where the input function is a unit impulse 
function /,(/). From the results of section 4.5 we have that the input 
spectrum F^co) = 1 . Hence for a network with a transfer function 
Z T (<o) we have for the output spectrum F 0 (co) 

F 0 (o)) = Z T (co) F^co) = Z T (co) 
and the output signal f 0 (t) 




Z T (co) exp (/co/) dcu 
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and so 

Z T (a>) = p f 0 (t) exp ( —jwt) dt 

J—oo 

Thus we have the important result that the transfer function of a 
network and its impulse response form a Fourier pair. 

The impulse response / 0 (/) of the network is the output when a unit 
impulse is applied at time /=0. If the impulse is applied at time /=r 
the output is/ 0 (r—r). Further if the impulse is not of unit magnitude 
but of magnitude fi(r) then the output is given by f\if) Any 

arbitrary input function can be made up of a series of impulses of the 
appropriate magnitude and the output due to an input function/ r (r) 
is the integral of the outputs due to inputs consisting of impulses of 
magnitude /i(r) occurring at time r. Hence the output function f 2 (t) 
is given by 

flit) = f°° /i(t)/o(/-t)^t=/ 1 (/)*/ 0 (0. 

J — 00 

This is the convolution integral. The output signal is therefore the 
convolution of the input signal with the impulse response of the net¬ 
work. The output spectrum is the product of the input spectrum and 
the transfer function of the network. We may compare these results 
with those obtained in section 3.13. 

5.8 The Method of Paired Echoes 

We shall now consider the effect on the output waveform of a filter 
which is not flat in the pass band. We have already seen that this will 
introduce distortion into the output and we now develop a method 
of calculating what this will be. In Fig. 5.4 we show the characteristic 
of an ideal filter by the dotted line. This has a value of one in the pass 
band and a sharp cut-off at co c = tt/t. We know that this characteristic 
must be an even function of co so the simplest form of variation we 
can introduce is a cosinusoidal variation with one complete cycle in 
the pass band. If this has amplitude a then the filter characteristic is 
represented by 

| Z T (co) | = 1 -f a cos cot — 7t/t < co < ttJt 
= 0 elsewhere 
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In order to avoid the effects of the sharp cut-off we suppose that 
the input signal has a spectrum whose components are negligibly 
small outside the pass band, i.e. for frequencies greater than 7 r/r. The 
ideal filter of Fig. 5.4 would then pass the signal without modification. 
Thus with this assumption any modification to the signal occurring 
when it is passed through the filter shown in Fig. 5.4 will be due 



A filter characteristic with 
amplitude variation in the 
pass band 

entirely to the ripple in the pass band. We shall also assume that the 
phase characteristic of the filter is ideal and from section 5.3 this is 
given by 

</>((x)) = — c oT 

Hence 

Z r (co) = (l4-tfCOScor)exp(— jwT) 

= exp ( —jcoT) 4- \a exp { -jw(T- r)} 4* \a exp { - ja)(T+ r)} 

Therefore if F}(cu) is the input spectrum, the output spectrum is 
given by 

W = F}(w)[exp(-ycur) + Jaexp{-yco(r-T)} 

+ ia exp { —j<o(T + r)}] 

We can now apply the shift theorem to each of these three terms 
and if /}(/) is the input signal then we have for the output signal 

foO) =fi(t-T)+Wi(t-T+T)+WtQ-T-r) 

This function is shown in Fig. 5.5. It is seen that in addition to the 
main signal which occurs after a delay T two smaller signals the same 
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shape as /■(/) but of amplitude \a occur spaced symmetrically about 
the main signal and at time r from it. These are termed paired echoes. 
In Fig. 5.5 these are shown separated from the main signal but in 
some cases they would overlap. The shape of the output signal would 
then be different from that of the input signal but it would still be 
possible to resolve it into the sum of three signals, a main one and a 
pair of echoes. Thus the distortion introduced by a cosinusoidal 
variation of filter characteristic can be represented by one symmetrical 
set of paired echoes. It is possible to represent any variation of filter 


fi(t) f 0 (t) 



Fig. 5.5 

Input and output of a filter with an 
imperfect amplitude characteristic 


characteristic in the pass band by means of a Fourier series. Each 
component will give rise to a pair of echoes and the resulting distortion 
can therefore be represented by a set of symmetrical paired echoes. 

We shall now consider the case of an imperfect phase characteristic. 
Since this must be an odd function of frequency it can be represented 
by a Fourier sine series and the simplest case is given by 

= —wT+bsintor 

Assuming an ideal amplitude characteristic we have 
Z T (co) = exp{-j(<*)T- b sin <dt)} 

Hence the output spectrum is given by 

F q ((d) = Fj(aj) exp { -j(coT- b sin cor)} 

Now from the theory of Bessel functions we have 

exp {y& sin*} = J 0 (b) 4- {exp 0’*) ~ ex P ( —/*)}/ 1 (b) 

4- {exp C/2*) 4- exp (-j2x)}J 2 {b) 4-... 
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where the /’s are the Bessel functions of the first kind. Suppose for 
simplicity that b is small so that the higher order Bessel functions can 
be neglected, i.e. we can write to a good approximation 

exp {jb sin tor} = J 0 (b) 4- {exp (Jojt) - exp (- jwr)}J x (b) 

Then we have for the output spectrum 

F 0 (oi) = F t (o>) exp ( -ju>T)[J 0 (b) + (exp (Jwr) - exp ( -j'ojt )}/, ( 6 )] 
= J 0 (b) exp ( —jo)T) +Ji(b) exp { -jo*(T- t)}F,(w) 

-J { (b) exp { —jco(T + t)} F t (w) 

Hence applying the shift theorem 

foU) = Ub)f i {t-T)+J l {b)f t {t-T+r)-J l {b)m-T-r) 

We can make the further approximation for small b,J 0 (b) = 1 and 
J,(b) = lb so that 

W) = m-T)+W,(t~ T+r)-lbfi(‘-T-r) 

The curves are shown in Fig. 5.6 and again we see that a pair of 
echoes appears, but in this case the second echo is reversed in ampli¬ 
tude. Hence a variation of the phase characteristic in the pass band 


ib 

o 

-ib 


fi(t) 


m 


A A A 




Fig. 5.6 

Input and output of a filter with an 
imperfect phase characteristic 


will give rise to distortion in the output waveform which can be 
represented by a set of anti-symmetrical paired echoes, the earlier one 
being positive and the later one negative. 
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CHAPTER VI 


Application to Wave Motion Analysis 

6.1 Introduction 

In the previous chapters we have been dealing with the Fourier pairs 
/(/) and F(cS) where the first is a function of time and the second a 
function of frequency. The use of Fourier analysis is by no means 
restricted to time-frequency pairs. If the functions f(x ) and F(y) 
form a Fourier pair then the product (xy) must be a dimensionless 
quantity. We shall begin this chapter by considering the case of 
diffraction in which the function/will be a space variable instead of a 
time variable. We then pass on to consider modulation. In all cases a 
close parallel will be observed between the phenomena arising in 
diffraction theory and modulation theory simply because the under¬ 
lying mathematics is the same. 

6.2 Diffraction 

Let us consider an absorbing screen AB as shown in Fig. 6.1 whose 
transmission coefficient at the point x is given by/(x). Suppose the 



Fig. 6.1 

Diffraction by a screen with 
variable absorption 

screen is illuminated with monochromatic radiation of wavelength A 
and consider the amplitude of the resulting wave in the direction 6 . 
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The contribution of the element dx at the point jc has amplitude 
proportional to f(x) and phase (27rsin0/A);c, or Kx where 

K = (27rsin0/A) 

Therefore if the incident radiation is represented by the complex 
quantity Rexp(jojt) the contribution in the direction 0 due to dx is 
given by Rf(x)expj(a)t- Kx). Hence the total contribution from the 
whole screen is given by 

J ^ Rf(x) expj(a)t-Kx) dx = J? exp CM) J* w f(x) exp ( - jKx) dx 

and the radiation in the direction 0, relative to the incident radiation 
may be written 

/(*) exp (-jKx) dx = F(K), K = 27rsin0/A 

Hence we see that/(*) and F(K) form a Fourier pair. Note that 
Kx = (27rsin0/A)x is dimensionless as required. The function f(x) 
represents the transmission characteristic of the screen. The complex 
function F(K) represents, in both magnitude and phase, the vector 
defining the radiation in the direction 0 (where 277-sin 0/A = AT). In 
this analysis we have assumed that no phase change is introduced by 
the screen and so f(x) is a real function. The function F(K) however 
does contain a phase term and so is complex. This case then corre¬ 
sponds to the most common case arising with the pair/(/) and F(co) 
in which/(0 is real and F(cu) in general complex. 


6.3 Diffraction at a Slit 

Consider the case of diffraction at a slit. Let the slit extend from 
x= 8 to at = +8 as shown in Fig. 6.2(a). Suppose the amplitude of 



f(4 

A 







-8 0 +8 


X -> 

(b) 


Fig. 6.2 

Diffraction by a single slit 
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the light transmitted by the slit is A and that the screen is completely 
opaque elsewhere. Then f(x) is as shown in Fig. 6.2(6). If we let 
AS = 1 then from the results of section 4.2 we can write down 
immediately 

sin KS sin (2778 sin 0/A) 

( “ K8 ~ (2tt8 sin 0/A) 

The intensity of the radiation is proportional to the square of the 
amplitude and hence varies like 

_ sin 2 (27r8sin0/A) 

(27rSsin0/A) 2 

which is a well-known formula in physical optics. 


6.4 The Diffraction Grating 


The formula for the diffraction grating can be developed in a similar 
manner to that for the slit. The case of the slit corresponds to a pulse 
and the grating corresponds to a finite train of pulses. Hence we have 
immediately from section 4.6 


/ = 



smKSsmNKd\ 
KS sin Kd ) 


2 


as the variation in intensity produced by a diffraction grating con¬ 
sisting of N slits of width 28 and spacing d. This again is a well-known 
formula in physical optics. 

It is of interest to note the effect of displacing the grating in the x 
direction. If the amount of the displacement is x 0 then f(x-x 0 ) 
represents the variation of the transmission characteristic with x. By 
the shift theorem this gives an amplitude variation exp ( - jKx 0 )F(K ), 
and the intensity is given by 

I exp ( -jKx 0 ) F(K) | 2 = |F(tf)| 2 

Thus the variation of intensity is not affected by displacement of 
the grating. 

We now consider the case of a sinusoidal grating. In this case the 
transmission characteristic varies sinusoidally across the grating and 
we have 


/( x) = 1 -f a cos (2ttx/1) 
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By the superposition theorem this will give rise to a transform 
consisting of the sum of the transform of 1 and of the transform 
of acos(27rx//). The first term corresponds to a beam in the direction 
0 = 0, i.e. an undeflected beam and the second gives rise to two 
components in the directions given by K = ± 2tt■//, or 

27rsin0/A = ±27r// 

Hence sin 0= ±A//. 

The beams in this case are infinitely narrow since we assumed a 
grating of infinite extent. In optical terminology the whole of the 
energy of the incident beam becomes concentrated in the zero and 
first-order spectra after transmission through a sinusoidal grating. 


6.5 Amplitude Modulation 

Suppose we have a waveform h{t) such that |/z(/)| ^ 1 for all t and let 
H(cu) be the Fourier transform of h(t). If now we modulate the 
amplitude of a carrier cos a) c t by h(t ) then we have the resulting 
waveform f(t) 

f(t) = {l+A(/)}costt> c f 
Expanding the function we have 

/(/) = cos a> c t+l exp (ja> c t) h(t) + i exp ( —jo) c t) h(i) 

Therefore the spectrum of the modulated carrier consists of two 
spectral lines representing the carrier plus, by the shift theorem, the 


w(«)t 


H(U >)t 


-<± 


<*)—> 


D- 


(Oc 


(*) (b) 

Fio. 6.3 

The development of the spectrum for amplitude 
modulation 


two terms <o + co c ) and \H (a> — a> c ). Suppose for simplicity H(w) 
is real, then as we have seen, it must be symmetrical about cu = 0 as 
shown in Fig. 6.3(a). The spectrum of the modulated wave is then 
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obtained by halving the spectrum and translating one half to the 
point a) = o) c and the other to w = - cu c . Note that from this point 
of view negative frequencies of H(w) after translation to the new 
origin will appear in the region of positive frequencies in the spectrum 
for the modulated wave. To complete the spectrum two lines must be 
added to represent the carrier. 

Suppose for example h(t) = cos pt. Then the spectrum is given by 
two impulses (or lines) of amplitude n (section 4.5) at frequencies 
co = ±p. The spectrum of a carrier when modulated by this function 
can be obtained by the construction just described. The result is 
shown in Fig. 6A(b) and the two familiar sidebands removed from 
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Fig. 6.4 

(a) The spectrum of a cosine wave and (b) a cosinusoidally modulated 

carrier 


the carrier by ±p appear. The spectrum F(w) of f(t) in this case is real 
and Fig. 6.4 can be considered as a plot of Re{F(co)} or of a{co) as 
defined in Chapter III. 

Now consider the case of sinusoidal modulation of a cosinusoidal 
carrier, i.e. 

h(t) = sinp/ 

and /(/) = (1 + sin pt) cos <o c t 

The complex spectrum H(a>) is now pure imaginary and 
Im {H(<o)} = -jb(<o) 

(section 3.4). We also know that the spectrum must be anti-sym¬ 
metrical. If we plot b(w) then we get two spectral lines, one of magni¬ 
tude + 7T at frequency +p and one of magnitude - 7 r at -p, as shown 
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in Fig. 6.5(a). Applying the construction developed above gives a 
spectrum shown in Fig. 6.5(6). Note that the construction leads to a 
spectrum which is also anti-symmetric as b(co) must always be. The 
complete spectrum for f(t) also contains a real part a(oj) consisting 
of two spectral lines of magnitude n and frequencies ± w c representing 
the carrier. 




(a) (b) 

Fio. 6.5 

(a) The spectrum of a sine wave and (6) a sinusoidally modulated 
carrier 


6.6 Amplitude Modulation and Paired Echoes 
In the last section we considered a time function which was modulated 
in amplitude sinusoidally. In section 5.8 we considered a function of 
frequency whose magnitude varied sinusoidally. We shall now 
demonstrate how these two effects are related by the principle of 
duality. 

Consider first the case of amplitude modulation of a sinusoidal 
carrier sine o c t. The modulated wave is represented by 

/(/) = {l 4- 6(0} sin o) c t 
= (i i 

{l + h(t)}j sin w c t = i{l + /i(/)}{exp (ja> c t) - exp ( -jw c r)} 

If 6(0 = cos pt then the spectrum H(wi) is real and hence the spec¬ 
trum of (1 +cospOysineo c r is also seen to be real. The spectrum can 
be obtained from the spectrum of cos pt and the shift theorem and is 
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shown plotted in Fig. 6.6(6). It may be noted that although this 
spectrum is real it is not symmetrical. The reason for this is that the 
symmetry of a(w) was deduced on the assumption that the original 
function of time was pure real and in this case this is not so. Figure 
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Development of response of 
filter with imperfect amplitude 
characteristic from amplitude 
modulation spectrum using the 
duality theorem 

6.6(a) shows the spectrum for the function (1 + cos/?/)cosa> c f, and 
the spectrum for the function 

/(/) = (l+cos/7/)(cosaj c /-ysinw c r) = (1 +cos/>/)exp(-ya> c /) 

is obtained by the superposition theorem and is shown in Fig. 6.6(c). 
Let this function be g(w). Then the signal of g( qj) is f(t), therefore by 
the principle of duality the signal of f(cS) is g(- /)/27 t. The expression 
for /(a>) is obtained from the expression for f(t) by replacing the 
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variable t by the variable a> and replacing the particular values of 
frequency p and o> c by the particular values of time r and T. Hence 

/(o>) = (1 + cos cut) exp (-ycuT) 

The signal of /(cu), i.e. g(- t)l2ir is shown in Fig. 6.6(d). The func¬ 
tion /(<u) represents the characteristic of a filter with an amplitude 
response (1 + cos cur) and a delay T. Hence the graph of Fig. 6.6(d) 
represents the output from such a filter for a unit impulse input 
(section 5.7), which shows that in addition to the main signal at time 
T a pair of echoes appear at T± r. (These results should be compared 
with those of section 5.8.) 

We see then that a cosinusoidal modulation of a time function gives 
rise to a pair of symmetrical sidebands, and a cosinusoidal variation 
with frequency of the amplitude response of a filter gives rise to a 
symmetrical pair of echoes. 

6.7 Phase Modulation 

Phase modulation is rather more difficult to deal with mathematically 
than amplitude modulation. It may be seen why this is so as follows. 
Suppose we wish to modulate the phase of a carrier of frequency cu c 
with the function m(t). Then the resultant phase of the wave is given 
by 

<P(t) = a> c t+m(t) 

and the signal is 

/(/) = cos {cu c t 4* m(t)} 

= i [ ex P Oc 0 exp {jm(t)} +exp ( - jeo c t) exp { -jm(t)}] 

The function m(t) now appears within the exponent and this makes 
the general case difficult to deal with. We shall therefore consider the 
case of sinusoidal modulation and for simplicity construct the case 
for f(t) symmetric so that the spectrum is pure real. Let 

f(t) = <f>{o) c t+m(t)} 

= 

r(t) = a> c t+m(t) 


where 
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We wish to construct the case for f(t) symmetric, i.e./(- /) =/(f). 
If f(t) is symmetric then <f>(r) is also symmetric 

<f>(7) = <K- r) 

Also, since /(/) = #r), /(-/) = <£(-r) 

Now /(-f) = <f>{-a) c t+m(-t)} 

and 4>(-t) = 

hence m(—t) = — m(t) and m(t) must be anti-symmetric. Thus since 
we wish to consider sinusoidal modulation we consider 

m(t) = <f) 0 sinpt 

The resulting modulated wave is then given by 
f(t) = cos {oj c t+</> 0 sin pt} 

= i[ ex P (M 0 ex P Oto sin pt) + exp ( - jo) c t ) exp ( -j<f> 0 sin pt)] 

Now from the theory of Bessel functions we have 

00 

exp(y<£ 0 sin/>r) = S J n (<f>o)^P(Jnpt) 

n— —oo 

exp(y^ 0 sinp/) = J 0 (<f> 0 ) + J , (<f> 0 ) exp (jpt) 

+J-i(<f>o) © X P (~jpO +• • • 

Wc can expand the expression exp(-y^ 0 sin/>/) by noting that 
exp {-y<£osin (/>/)} = exp {j<f > 0 sin (-pi)} 

exp{—y^ 0 sinp/} = J 0 (4> 0 ) + J i (4>o) exp( -jpt) 

+/_ i (<£o) exp (jpt)+... 

Hence expanding the expression for/(f) we have 
fit) = i[/o(^o) exp (jco c t)+J i(^ 0 ) exp {j(<o c +/»)<} 

+J- i(to) exp {j(o> c -p) /}+... 

+ Mto) exp ( -j<» c f)+/, (<t> 0 ) exp {-j(to c +p) /} 

+/-i (<f> 0 ) exp { -j(w c —/») /}+...] 

'•AO = /o(^o) cos f+/, (<f> 0 ) cos(w c +p)t+J_ l (4> 0 ) cos (o} c -p) t 

+ M’t’o) cos (<o c + 2p) t +/_ 2 (<f> o) cos (a» c - 2 p)t+... 
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This expression may also be written 

00 

/(') = 2 Uto) C0S ( w c + n P) t 

n— —oo 

Thus we see that for sinusoidal modulation we get an infinite series 
of spectra spaced at multiples of the modulation frequency above and 
below the carrier. If there were more than one modulating frequency 
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Fig. 6.7 

The spectrum of a phase modu¬ 
lated carrier (diagrammatic) 


then it may be shown that in addition to the two infinite sets of 
sidebands corresponding to the two frequencies, all the beat frequen¬ 
cies also occur. This is why the general case is so complex. 
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Fig. 6.8 

The spectrum for small phase 
deviation 

If the maximum phase excursion <j> 0 is small then we can sketch the 
spectral lines bearing in mind that /_„(*) = (- l) n J„(x). The sketch 
is shown in Fig. 6.7. The way in which the spectrum develops as <f> 0 
becomes larger is discussed in works on phase modulation and 
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frequency modulation (see, for example, Arguimbau and Stuart 
(1956)). lf<t> 0 is sufficiently small then / 0 (^o) = 1 and/j^o) = i^o anc ^ 
we may neglect all higher order terms. The expression for fit) then 
becomes 

fit) = cos <D C t + o cos (a> c +/>) t - \<f>Q cos (w c -p) t 

This expression corresponds closely to that for amplitude modu¬ 
lation, the difference being that the lower sideband is negative. The 
spectrum for this case is shown in Fig. 6.8. 

6.8 Phase Modulation and Paired Echoes 

In the same way as we related amplitude modulation to distortion of 
a filter amplitude characteristic so we can relate phase modulation to 
distortion of a filter phase characteristic by the principle of duality. 
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Development of response of 
filter with imperfect phase 
characteristic from phase modu¬ 
lation spectrum using the 
duality theorem 

Consider the case just discussed in which <f> 0 is small so that only 
terms of the first order need be taken into account. Following the 
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method used for the amplitude case we form a spectrum in the 
negative region of co. Thus we consider 

/(/) = exp{-yco c /+y^ 0 sin/?r} 

= M<f>o) exp ( -jw r t)+J l (<f> 0 ) exp ( -jw c t +jpt) 

+ J - exp (-]<*> c t-jpt) 

= exp ( -ju c t) + i<f> 0 exp { -j(w c -p) /} 

- Ho exp{-j(o> c +P)*} 

Hence the spectrum g(cj) is as shown in Fig. 6.9 (a). By the principle 
of duality if/(/) is the signal of g(w) then the signal of/(a>) is 

g(—t)/27T 

Hence the signal shown in Fig. 6.9(6) has spectrum 
/(co) = exp { -j(wT- <f> 0 sin cor)} 

Hence Fig. 6.9(6) is the output of a filter with this characteristic as 
transfer function (i.e. a distorted phase characteristic) when the input 
is an impulse. Again, in addition to the main impulse at Ttwo echoes 
occur, the later one being inverted. (These results should be compared 
to those obtained in section 5.8.) 

6.9 Vector Representation 

In the preceding chapters we have used various methods of represent¬ 
ing the functions we have been discussing. We have used analytic 
expressions to represent functions of time directly; we have trans¬ 
formed these into functions of frequency, which in general were 
complex, using either an analytic expression or alternatively plotting 
the real and imaginary parts of the spectrum graphically; and we have 
used the idea of vectors (which vary with time) plotted on an Argand 
diagram. In this section we shall consider the vector representation 
of modulated waves. 

We saw in section 3.6 that a cosine function can be considered as the 
sum of two rotating vectors representing the functions £exp(ya> c /) 
and £exp(— jo) c t). These vectors produce a real variation cosu> c /. 
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The relative phase of the vectors is obtained most readily by con¬ 
sidering t = 0 and at this instant they both lie along the positive real 
axis. The spectrum is pure real and consists of two lines at ± w c of 
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(a) Argand diagram for vector representation and (6) and (c) the 
spectrum of a cosine wave 


magnitude rr. This can be summarized diagrammatically as shown in 
Fig. 6.10. Figure 6.10(a) is the Argand diagram drawn for t = 0. The 
two vectors should be coincident but they are drawn slightly separated 
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for clarity. Figures 6.10(6) and 6.10(c) show the plots of a(oj) and 
6(<o), the latter being zero everywhere in this case. 

The case for sin u) c t is shown in Fig. 6.11. In order to produce a real 
variation of sin w c t f i.e. a variation along the x axis, we multiply the 
positively rotating vector by —j so that at time / = 0 it lies along the 
negative y axis and the negatively rotating vector by +/ At time t — 0 
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the two vectors clearly cancel and the resultant is zero and in general 
we ha ve £/{exp ( -jw c t)- exp (jw c /)} or {exp (jw c t) - exp ( - jw c t)}/2j, 
which gives a real variation sinco c r as required. The spectrum is 
therefore pure imaginary and is plotted in Fig. 6.11® and (c). 

We shall now develop a series of these diagrams to illustrate various 
cases of amplitude and phase modulation. We begin with the simplest 
case of cosinusoidal amplitude modulation of a cosinusoidal carrier. 

fit) = (1 + a cos pt) cos w c t 

In this case the carrier is represented by a pair of vectors rotating at 
± (d c as before and in addition each has a pair of sidebands. These can 
be considered as rotating at rates ±p with respect to the carrier vector. 
Thus the sideband vectors can be drawn from the tip of the carrier 
vector and must then be considered to be carried round with it. At 
t = 0 the situation is as shown in the Figure. ( Note: The vectors are 
actually coincident with the real axis but are shown separated for 
clarity.) The sidebands at this instant combine in such a way as to 
increase the magnitude of the resultant vector by the maximum 
amount, as required by the analytic expression. At all times the 
sideband vectors combine in such a way as to produce a resultant in 
the same direction as the carrier vector they are associated with and 
hence their effect is to modify its length or to modulate its amplitude. 
The resultant carrier vectors are always identical in length and 
combine to produce an overall resultant along the real axis. We can 
also consider this construction as consisting of producing first a 
cosinusoidal carrier by means of the two vectors rotating at ± w c and 
then the cosinusoidal modulation of their amplitude at frequency p 
by constructing a term acospt with the carrier vector regarded as real 
axis. Thus at / = 0 the vectors rotating at ±p lie along the carrier 
vectors, for this case along the real axis of the Argand diagram. Since 
at / = 0 all vectors lie along the positive real axis the spectrum is pure 
real and all terms are positive. The plots of a(o>) and b(w) are shown 
in Fig. 6.12. 

The second case is when the modulation is sinusoidal and the 
resulting analytic function becomes 

fit) = (1+ flsin/?0cosc*> c / 
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This case is shown in Fig. 6.13. The same considerations apply in 
this case. The vectors representing the sidebands are constructed at 
the tip of the carrier vectors and at / = 0 combine in such a way as to 
leave the length of the carrier vector unmodified. Again the sideband 
vectors always combine to produce a resultant in the same direction 
as their associated carrier and the two carriers combine to produce an 
overall resultant along the real axis. Alternatively, we can regard this 
construction as before as first producing the carrier by means of two 
vectors rotating at speeds ±co c and then producing sinusoidal 
modulation of their amplitude by constructing the asinpt term with 
the carrier vector regarded as real axis. Thus at / = 0 the sideband 
vectors all lie parallel to the imaginary axis in the Argand diagram. 
Hence the complete complex spectrum has a real part consisting of 
two lines at ± cj c and an imaginary part consisting of four lines at 
( w c± p) and (— io c ± p). These are shown in Fig. 6.13(6) and (c). 
Note that the sideband vector in Fig. 6.13 (a) associated with the 
+ o) c vector and rotating positively, i.e. corresponding to the sideband 
(cn c + p) is drawn in the negative direction and hence corresponds to a 
positive line in 6(o>) since 6(co) = — Im {F(o>)}. We can also obtain 
the spectrum from the construction developed in section 6.5. The 
modulating waveform is sin pt and its spectrum consists of two lines 
as shown in Fig. 6.11(c). We halve the magnitude of these lines and 
translate the result to the new origin which gives the same result as 
before, shown in Fig. 6.13(c). To complete the spectrum the real part 
consisting of two lines representing cos co c t must be added. 

Two more cases arise for amplitude modulation, when the carrier 
sinusoidal. The analytic expressions are (1 +flCosp/)sina> c f and 
(1+Asinp/)sinai c /. These may be treated in a similar way and the 
results are shown in Figs. 6.14 and 6.15. 

We can consider the four similar cases for phase modulation. For 
simplicity we shall suppose that the maximum phase modulation <f> 0 
is small so that only the first pair of sidebands need be considered 
(see section 6.7). Considering first the case f(t) = cos(o) c /+<£ 0 sin/?/) 
we must produce a variation at right angles to the carrier vector 
in order to modulate its phase. Hence the sideband vectors are as 
shown in Fig. 6.1 6(a) and it is seen that at t = 0 all vectors lie along 
the real axis and the spectrum is therefore pure real. This shows 
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diagrammatically why in order to produce a symmetrical spectrum 
(i.e. a pure real spectrum) in the case of phase modulation, the modu¬ 
lating waveform must be a sine function as we proved analytically in 
section 6.7. A small value of <f> 0 means that the ratio of the length of 
the sideband vector to that of the carrier vector is small. In this case it 
can be seen that the effect is to modify the phase of the vector without 
affecting its amplitude to a first approximation. As <f> 0 becomes 
larger, however, the effect on the amplitude becomes more important 
and it can be seen that the resultant vector is larger at the maximum 
phase excursion in both positive and negative directions. Hence an 
amplitude modulation is introduced at twice the frequency. This can 
be cancelled out by a suitable term at frequencies co c ±2p and this is 
how the second-order sidebands arise. These in turn modify the phase 
as well as the amplitude of the resultant calling for higher order 
sidebands to correct the effect. Thus we can see from the vector 
representation why a simple sinusoidal modulation of phase gives 
rise to an infinite series of sidebands. 

Three more cases arise for phase modulation and the vector dia¬ 
grams corresponding to the analytical expressions cos (co c t + <f> 0 cos pt ), 
sin(co c /+<£ 0 sin/7/) and sin(a7 c /+<£ 0 cos/?/) are shown in Figs. 6.17, 
6.18 and 6.19. The spectra are readily obtained from the vector 
diagrams remembering that b(ui) = -Im{F(ci7)}. 

6.10 Modulated Pulse Trains and Grating Faults 
We have already considered the spectrum of a train of pulses 
(sections 2.6 and 4.6) and in the following sections we shall investigate 
how this spectrum is modified when the pulse train is modulated. We 
can modulate a pulse train in various ways, and we shall consider 
pulse amplitude modulation (P.A.M.), pulse phase (or position) 
modulation (P.P.M.) and pulse width modulation (P.W.M.). In the 
first two cases the simplest results are obtained if a train of impulses 
is used. Such a simplification is not possible however in the case of 
P.W.M. When pulses of finite width are introduced it is necessary to 
define the modulation process with some care, as will be shown. 
Using the results obtained at the beginning of this chapter it can be 
shown that the spectra of modulated pulse trains will be the same as 
the optical spectra produced by a diffraction grating with various 
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types of periodic faults in the ruling. The two sets of phenomena will 
be related in subsequent sections. 


6.11 Pulse Amplitude Modulation 

We shall consider a train of impulses modulated in amplitude by a 
cosine function as shown in Fig. 6.20 (a). Suppose the unmodulated 



Fig. 6.20 

(a) The pulse amplitude modulation signal and (b) its spectrum 


pulse train is /(/) and has spectrum F(o>). Then the modulated train 
is given by 

h(t) = (1+ a cos pt)fit) 

= fit) + ha exp 00/(0 + ha exp (-00/(0 

Hence by the shift theorem the spectrum of the modulated train 
is given by 

H{w) = F(a>) -f h a F(cu ~P) + h a F(w -f p) 

The resulting spectrum is shown in Fig. 6.20(6). The value of 
cuJItt for the magnitude of the spectra is that appropriate to a train 
of unit impulses occurring at intervals 2tt /co r . In this case of course 
the height of the pulses would be infinite. If on the other hand we 
consider impulses of finite height then the magnitude of the spectral 
lines is zero as we saw in section 2.6. To produce a case in which both 
the time and frequency functions are finite we must consider pulses of 
finite width. In this case the spectral lines are modified in amplitude 
by the term sin(fca> r T)/(fcT) where k is the order of the spectrum and 
2 t is the pulse width. The function h(i) was obtained by multiplying 
fit) by the factor (1 + acos/?/) so that at each instant of time/(/) is 
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modified by the value of (1 +acos/?/) at that same instant. Thus the 
resulting pulses are not rectangular pulses and the modulation 
modifies their shape as well as their amplitude. 

From the spectrum shown in Fig. 6.20(6) we see that a low pass 
filter can be used to separate the term at o> = p from the rest of the 
terms and hence produce the original modulating wave. Thus a low 
pass filter acts as a detector. The spectrum also shows that the nearest 
term from which the modulation term is to be separated occurs at 
(w r - p). Hence for separation to be possible we require 

<*> r -P > P 

or a) r > 2p 

The pulse repetition frequency must be at least twice the highest 
modulation frequency to be transmitted, otherwise the lower side¬ 
bands of the first harmonic of the repetition frequency overlap with 
terms occurring at the modulation frequencies and it becomes 
impossible to recover the original modulating waveform without 
distortion. 

6.12 Ghost Spectra 

The function h of the last section could represent a diffraction grating 
in which the absorption of the transparent portion had a sinusoidal 
variation with distance. Such an effect is produced by a variation in 
the heaviness of the ruling of the lines. The function H would then 
represent the resulting spectrum produced, as we saw in section 6.2. 
Each spectral line is accompanied by two fainter lines equally spaced 
on each side of it. These have been termed ghosts. The intensity of the 
ghost produced by this type of fault is independent of the order of 
the spectral line it accompanies as can be seen from the plot of H 
in Fig. 6.20(6). 

6.13 Pulse Phase Modulation 

We shall first consider the simpler case of a train of impulses. We may 
write the expression for these 

oo 

/(') = ^ ^ eX P U k<t >r>) 

k *= —oo 
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We now, by analogy with the method used in the derivation of 
phase modulation, replace (o r t by (< w r t+<f> Q sinpt ). Hence we have 

00 

hit) = ^ 2 exp *+&> sinp/)} 

k=* -oo 
oo 

O) v - ' 

= 2^ CX P 0 ko >r 0 eX P OMo sin />0 

k=- 00 


and using the relation from Bessel function theory 

00 

exp(jk<f> 0 sinpt) = £ J„(k4> 0 )exp(jnpt) 

n=—ao 

we have 

00 oo 

AO = ^ 2 eX PO ko> rO 2 / «(*^o)«*P(/'»PO 

k=— oo * /!=■— oo 


= ^ 2 2 exp {X^r+«/>)'} 

k =-oo n = —oo 

We can simplify this as follows. If k = 0 then J n {k <!> 0 ) = J n (0) = 0, 
n ^ 0. Hence for k = 0 the second summation consists of only one 
term, viz. / 0 (0)exp(0) = 1. We can therefore write 


.no =^[i+ 2 (2 + 2) Wo) ex p ix km r + np) r}l 

[ CO 00 

i+ 2 2w^ exp ^ + ^' 

n = — oo k = l 

+J„(-k<l> 0 )expj(-km r +np) /}J 

[ 00 00 

1+2 2 eX PX* w r + n P) t 

n= —oo k—l 

+J-„(-k<l> 0 )exp[-j(ka> r +np)t]}^ 
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and since 

J-ni-x) = J n (x) 

[ 00 00 
i+ ^ 

n— —oo *=1 

+ exp [ -j{kw r + np) /]}j 

[ 00 00 

1+2 ^ ^ C0S (^ w r + »/>) ‘ I 
Expanding the summation terms we have 

/CO = -^[i+Jo(<l > o)cosa>rt+J l (<l> 0 )cos(a) r +p)t 

+J -1 (^o) cos ( w r “/*) t + • • • 

+/o(2^o) 008 2o> r t +/, (2<£ 0 ) cos (2a> r +/>) / 

+/_ i(2^ 0 ) cos (2o> r -/>) 

Figure 6.21 is a diagrammatic sketch of the spectra for small 
deviation <f> 0 . It consists of a series of lines at the harmonics of the 
repetition frequency each having an infinite set of sidebands similar 
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The pulse phase modulation spectrum (diagrammatic) 

to those for a phase modulated sinusoid. Unlike pulse amplitude 
modulation the sets of spectra occurring at w n 2w n etc. are not 
similar. The magnitude of the sidebands forming the set centred at 
kw r being determined by the set of functions J n (k<f> 0 ). For small <f > 0 , 
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the sidebands are larger and more of them will be in evidence the 
higher the harmonic of the repetition frequency, as shown in Fig. 6.21. 
As </>o becomes larger the spectra will develop in a manner analogous 
to that for a phase modulated sinusoid. In the case of pulse phase 
modulation considered in this section there is no spectral term at the 
modulation frequency p. Hence again in contrast to pulse amplitude 
modulation a low pass filter will not act as a detector. 

6.14 Rowland’s Ghosts 

The results of the last section can again be interpreted in terms of the 
diffraction grating. In this case the function f(x) would represent a 
grating with a periodic displacement of the lines from their true 
position. Such faults occur as a result of slight imperfections in the 
ruling engines which can take various forms. An imperfection in the 
screw gives a periodic error with a period equal to the pitch of the 
screw. In this case the ghosts are known as Rowland’s ghosts. These 
are more troublesome the higher the order of the spectrum as may be 
seen from Fig. 6.21 and in the high orders sometimes results in 
fuzziness of the spectral line. If <f> 0 is small then we can expand 
J„(k<f > 0 ) as a power series and neglect all but the first term, so that 
Wo) = Wo)72"«!. The amplitude of the ghost nearest to the 
main line will therefore be proportional to k, and the intensity to 
k 2 , i.e. the intensity of the ghost varies as the square of the order of 
the spectrum. 

6.15 Pulse Phase Modulation—Further Consideration 

In section 6.13 we obtained the formula for a phase-modulated train 
of impulses by replacing a) r t by (a) r t+</> 0 sinpt). Thus if the train of 
impulses is represented by /(Owe obtained the function/(f+ t 0 smpt) 
where t 0 = ^ 0 /w r . Suppose the impulses occur with a spacing in time 
of T. Then/(/) is zero except at the points t = f, where = ±nT and 
n = 0, 1,2..., and at these points an impulse occurs. If now h(t) is 
given by/(/ +/ 0 sin pi) then h(t) is zero everywhere except at the points 
t t where / / +r 0 sinp/ / = ±wTand at these points an impulse occurs. 
Hence it can be seen that the value of the modulating waveform 
affecting the position of the pulse is its value at time //. Thus the 
position of each pulse is determined by the value of the modulating 
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function at the instant at which that pulse occurs. The spectrum 
derived in section 6.13 corresponds to a pulse train modulated in this 
manner. 

However the process used in section 6.13 is not justified because 
the function h(t) representing the phase-modulated train is not in 
general periodic. It will only be so if w r /p is rational. Hence in general 
we are not justified in assuming a Fourier series type of representation 
as we did in section 6.13 and we now consider a derivation in which 
this assumption is not made. In section 4.6 it was shown that a scries 
of N pulses of unit area and width 2r occurring symmetrically about 
the origin could be represented by the Fourier integral 
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jl r 

"27rJ_ a 


sin cut sin \No)T 


exp (jwt) dco 


t sin 

We now obtain h(t) by replacing / by t+t 0 smpt to give a phase- 
modulated pulse train of finite duration 

KO =fO+t 0 sin/?/) 


jl r si 


sin cut sin \Nu>T 


exp {ja)(t+1 0 sin pt)} daj 


2 7rJ_oo cut sin i<oT 

Examining the limit N-* co by the same procedure as we used in 
section 4.6 we see that, in the limit, the only contributions to the 
integral will come from the points where sin icuT= 0, or the points 
cu = kw r where cu r = 2tt/T. Hence we can write 
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exp {jka> r (t + 1 0 sin/?/)} 


Now letting t->0 to produce impulses we obtain 
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6.16 Phase Modulated Pulses of Constant Width 


We saw in the last section that the function obtained by replacing 
t by (f + / 0 sin/?/) produced a series of impulses whose time of occur¬ 
rence was affected by the value of the modulating function at the 
instant when the pulse occurred. Similarly if the original function had 
represented a series of steps instead of impulses the time of occur¬ 
rence of the step would be affected by the value of the modulating 
function at the instant the step occurred. Now if the pulses are of 
finite width we can consider that they are composed of two steps, the 
first a positive one and the second a negative one occurring at a time 
2t, the width of the pulse, later. Hence from the argument just given 
if we replace t by (/+f 0 sinpr) in the expression for pulses of finite 
width then the position of each edge of the pulses is modified by the 
value of the modulating function at the instant when the edge occurs. 
Since in general the modulating waveform will not have the same 
value at the trailing edge of the pulse as it had at the leading edge of 
the same pulse, such a pulse train will be width-modulated as well as 
phase-modulated. If then we require the spectrum of a train of phase- 
modulated pulses of constant width we must use a process which 
treats the leading and trailing pulse edges independently. 

In section 4.8 we showed that a unit step could be represented by 
the expression 




expO’^O 

jw 


dto 


Hence a series of AT unit steps occurring at / = - N'T 9 — ( N ' -1 )T ,... 
-T, O, T,... (N'- l)r, N’T, where N= (2 N'+ 1), will be given by 
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We now phase modulate the position of these steps by replacing 
/ by (/+ t 0 sinpt) giving 
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It was shown in section 4.6 that 


iV 

2 


exp (JncoT) = 


—N* 


sin \N(dT 
sin Jc oT 9 


N= 2AT+1 


Hence the expression above becomes 

iAr+ i_ p 0 exp +1 0 sin/?/)} sin jNcoT ^ 

ZnJ- oo y<o sinicoT ^ 

To produce a series of phase-modulated pulses of constant width 
2r we shift this waveform r to the left by writing (/+r) in the expression 
instead of t and subtract from it the same expression shifted t to the 
right. Hence we obtain 


1 f°° 1 
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Taking the limit N-+oo as before we get contributions to the integral 
only at points where sin^coT = 0 or co = k(o r where co r = 27r/Tand 

m2 2 w« r (k “- k +np)T wi(ic<>,+*P)’ 

k= —oo n= —oo 

which is the expression required in complex form. We can simplify 
this expression as we did the one in section 6.13. For k = 0 we have 
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For the rest of the series we have 
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S [w°) S1 ° ~ - {exp j(k< 0 , + np) t 

ft — —CO *=l L K 

+expy( - km r - np) /}j 

“ n 2-, 2 Wo)- r -cos (kw r + np) i 

n— —oo *=I K 

Hence for the complete expression we have 
W) ~ — |^> r t+ sinpr cos pt 

n=-oo k= 1 K J 

We see that there is now a term at the modulation frequency, p, 
i.e. ^osinprcosp/. However its magnitude, <f> 0 smpr, is a function of 
the modulation frequency and hence for general modulating wave¬ 
forms the low frequency terms do not represent the original modulat¬ 
ing waveform. A low pass filter would not therefore be a suitable 
detector. 

00 

The expression 2 Jft(k<f> Q )cos(ka) r +np)t represents a phase 

fl= —00 

modulation of the carrier cos ko r t by the function cos pt with a phase 
excursion of k<f> 0 (see section 6.7). A phase detector centred at k<o r 
would produce the original waveform from this set of terms. How¬ 
ever in the expression above this term is modified by the factor 
s\n(ka> r +np)T which again is a function of /?, the modulation 
frequency, so again we cannot obtain an undistorted output. 

It is common practice to detect pulse phase modulation by first 
converting it into pulse width modulation. 

6.17 Pulse Width Modulation 

In pulse width modulation the time duration of the pulses is modified 
by the modulating waveform. We shall assume the leading edges of 
the pulses are independent of the modulation and therefore occur at 
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constant time intervals and that the position of the trailing edge is 
modified in accordance with the modulation waveform. It can be 
seen that on this assumption it is not possible to construct a sym¬ 
metrical waveform. We shall suppose that a positive going edge (i.e. 
a leading edge) occurs at t = 0. Constructing the pulses by a process 
similar to that used above we have for the series of steps representing 
the leading edges 

J_ p expCAoQsin \N»T 

2nj J _ 00 o siniwT 


If the mean length of the pulse is r then the negative steps must first 
be displaced in time by +r which is obtained by replacing / by 
(/—r). We then wish to modulate their position (or phase) by r 0 sin/?/ 
which we do by replacing t by (/+r 0 sin/?/). The maximum pulse 
width will be (t+t 0 ), and the minimum (t-t 0 ). The depth of 
modulation will be r 0 /r. Carrying out the two operations above we 
have for the required series of negative steps 

wt 1 f 00 ex P “ T + T o s » n P*)} sin i Na>T r 

2 ■‘Y ~ I ;—■—~~~ (It*) 

2ttj J.qo c*> sm io)T 

Hence for the complete waveform we have 


1 f 00 1 

h(l) = — I - [exp (/of) - exp [jw(t - 


+ t 0 sin/>/)}] 


sin iNojT 

- do 

sin \oT 


taking the limit N^oo as before we have 



^ f e*p C/froy/) 

k— —oo ^ 


exp (jka> r t) exp {—jkco r r) exp (jkw r r 0 sin pt)) 

k J 

1 r exp(y^co r l) 

2tt/ Z [ k 

k= —oo U 


exp(jkco r t) 


00 1 

exp(— jkw r r) 2 *^«(^ <u r , ‘o) exp (Jnpt) I 
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which is the required expression in complex form. A sketch of the 
waveform represented by this expression is shown in Fig. 6.22. 
Simplifying as above we consider first k = 0 giving 

1 [exp (jkw r t) exp (jka» r t) , xfr/ ^ 

2nj\ - k - k - ex P ( ~jkui r r){7 0 (0) 

+ Ji(kw r r 0 ) exp ( jpt ) + /_ i(koj r r 0 ) exp ( -jpt )}| 


ht 


\k=0 


1 [exp (jkw r t) exp (jkw t) 

- i -*— 

exp(jktu.t) . e 

-£-exp ( —jk(o r r){ikot r r 0 exp (Jpt) 

— $ko> r r 0 exp(—jpt)} 

Jk=0 

1 [ exp {jk<x> r (t — Jt)} f 

2nj [ - k -l exp 0± koj r T ) - exp ( -j\kw r r)} 

- fa, t 0 exp {jkojjt - r)} (exp (Jpt) - exp ( -jpt)} 

J*=o 

1 [sin (\ko) r r) ] 1 

= ^-Jco r T 0 sinp/| = —\to r T-u) r T 0 sinpt] 
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Pulse width modulated signal 
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For the rest of the series (k ^ 0) we evaluate the two terms 
separately. The first gives 


L/V + y \ exp(jkco r t) _ 1 ^ exp (jkio r t) ex p (-jko> r t) 

nJ \4i k l7ri & 


2ir h ctr k 


1 -XT' sin(A:<u r f) 

n li k 

k— 1 


and the second gives 


ij{2* 2 | 

k— 1 k * 1 n = ■" co 

-ir'i i 


k= 1 n= — oo 1 


_ exp{ jktoXl t )) T o j exp )J 

-is2 2 h><y->W ,)c»P(^> 

7 1 «=-00*- 
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- 1 £ - J^ n (-kw r T 0 )cxp(-jnpt) 

— i 00 00 n 
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Hence we have finally 


1 r °° 1 

KO = — |a> r ( T - T oSinpf) + 2 ^ - sin (km r t) 


00 00 1 

“ 2 2 2 k Jn ^ kWr T °^ sin !(**“' +np ) t ~ ka} ' T 

k= 1 n=— oo 

This is a more awkward looking expression than those we have 
obtained previously because the function is not symmetrical about 
the origin. It is worth while making a check to see that it reduces to the 
correct form in the absence of modulation, i.e. r 0 = 0. In this case we 
have 


KO = “|oj r T +2 ^ sin (kco r /) — 2 ^ ^sin£a> r (/-T)j 

1 T °° 1 1 

= 2^U f r+2 2^{sin(Arw,./)-sinA:a> r (/-T)} 

L k— 1 J 

i r 00 i i 

= —1 0J r T+4 ^ - sin (Ik to, t) cos ku> r (/- }r) 

L k= 1 J 

We can check this from the formula developed in section 2.6. From 
these results we see that for pulses of duration r (note that in section 
2.6 a duration of 2r was used) occurring symmetrically about the 
origin we have 


fl o=r = 


2 0) r T 
2n 


(i) r T 
7r 


a k = — sin(JA:a > r r) 


and 
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so we can write 


/(f) = ia 0 + 2 a k coskw r t 

k=l 

CO T 00 2 

= ito + 2 dfe 8 ™ 1 ^ ku>rt ) cos ( kojf 

k=l 

ip ® i 

= —I w r T+4 ^ - sin (jigo r r) cos (kw r /) 

L /c = 1 J 

Now to obtain the same function as above we must shift this pulse 
train so that a leading edge occurs at t = 0 instead of t = - Jr. Hence 
replacing t by (f - ir) we have 

i r 00 

/(f) = ^|^r T + 4 ^ ^ SinT ) 008 WJ 

which checks the result above. It can also be seen in a general way 
from Fig. 6.22 that the component sin/?/ will be negative since the 
coefficient will be given by the formula 

2 f r/2 

~ /(f) sin/?/ dt 

1 J-772 

and it can be seen that the wide pulses occur where the sine function 
is negative and hence the negative contributions to the integral will 
be the larger. This is also in agreement with the expression obtained. 

In the expression for the pulse width modulated wave there is a term 
at the modulation frequency sin/?/ whose magnitude is independent 
of p and proportional to r 0 the depth of modulation. The undistorted 
modulation waveform can therefore be obtained from this and a low 
pass filter will act as a detector. 

6.18 Abbe Theory of Image Formation 

In the early part of the chapter we developed the idea of Fourier pairs 
in which the functions were functions of space variables instead of 
time variables and the similarity between diffraction patterns and the 
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frequency spectra of time functions was demonstrated. The Abbe 
theory of image formation is closely parallel to the case of filtering 
which we discussed in section 5.2. There we started with a function 
of time which we transformed into its spectrum, the spectrum was 
modified by the filter and the output time function was the inverse 
transform of the resulting frequency function. The Abbe theory of 
image formation states that for a truthful image of a structure to be 
formed by a lens the aperture must be wide enough to transmit the 
whole of the diffraction pattern produced by the structure; if only a 
portion is transmitted then the image produced corresponds to an 
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The formation of an image by a lens 


object whose complete diffraction pattern consists of the portion 
transmitted. Consider a one-dimensional structure represented by 
/(*). Then as we have seen the diffraction pattern is F(K) where 
K = (27rsin0)/A. If the aperture of the lens is limited then light rays 
within only a finite range of 6 are transmitted and hence F(K) is 
modified to F X (K) say, and the image becomes /,(*) the inverse 
transform. 

Consider for example a grating in the plane OB (Fig. 6.23) illu¬ 
minated by monochromatic light. Bundles of rays emanating from 
three of the slits are shown in the figure. The rays marked with one 
arrow are parallel to the axis and will therefore come to a focus at 
the point F { the principal focus of the lens. Suppose the rays marked 
by a double arrow are in the direction of the first spectrum. These are 
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all parallel to each other and therefore come to a focus in the focal 
plane FF' at some point F 2 . Similarly for the triple arrowed rays in 
the direction of the first spectrum on the other side of the axis. Thus 
the spectra are formed in the focal plane of the lens. The bundles of 
rays from the points in the object plane come to a focus in the image 
plane IM and form an image. If however a stop were placed in the 
focal plane so that only the first order spectra were transmitted then 
the image would not be a true image of the grating in the OB plane but 
would correspond to the inverse transform of the spectra transmitted. 
In this case where only the first order are transmitted a sinusoidal 
variation of intensity would be formed in the IM plane. The effects 
of passing only a limited number of spectra can be demonstrated 
experimentally. These results are of practical significance in micro¬ 
scopy where on occasions the object may differ considerably from 
the image due to modification of the spectra by the aperture of the 
instrument so that care is needed in interpreting the results. 
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APPENDIX 


DIFFERENTIATION OF A DEFINITE INTEGRAL 
WITH RESPECT TO ITS UPPER LIMIT 

We wish to evaluate 

d C at 

d,j„ Mdx 

where a is a constant. 

Let h(x) = jg(x)dx 

so that 

^A(at) = s(jc) 

Then 

j* g(x)dx = h(at)-h( 0) 

d C al d 

/.-Jo g(x)dx = -Kat) 

d r°' 

"Jtjo g w dx = ag ( af > 
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